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ABSTRACT 

We describe the Multi- Epoch Nearby Chister Survey (MENeaCS), designed to measure the cluster 
Type la supernova (SN la) rate in a sample of 57 X-ray selected galaxy clusters, with redshifts of 
0.05 < z < 0.15. Utilizing our real time analysis pipeline, we spectroscopically confirmed twenty-three 
cluster SN la, four of which were intracluster events. Using our deep CFHT/Megacam imaging, we 
measured total stellar luminosities in each of our galaxy clusters, and we performed detailed supernova 
detection efficiency simulations. Bringing these ingredients together, we measure an overall cluster 
SN la rate within i?2oo (1 Mpc) oi M2tofwtom8 SNuM (o.049+"°^<'+°™^ 



0.014-0.004 SNuM) and a SN la rate 



within red sequence galaxies oi OMltomltomo SNuM (0.04llg;°}^tg ™| SNuM). The red sequence 
SN la rate is consistent with published rates in early type/elliptical galaxies in the 'field'. Using our 
red sequence SN la rate, and other cluster SNe measurements in early type galaxies up to z ^ 1, 
we derive the late time (> 2 Gyr) delay time distribution (DTD) of SN la assuming a cluster early 
type galaxy star formation epoch of Zf = 3. Assuming a power law form for the DTD, \I'(t) cx , 
we find s — — 1.62 ± 0.54. This result is consistent with predictions for the double degenerate SN la 
progenitor scenario (s ~ ~1): and is also in line with recent calculations for the double detonation 
explosion mechanism (s ~ —2). The most recent calculations of the single degenerate scenario delay 
time distribution predicts an order of magnitude drop off in SN la rate ^^6-7 Gyr after stellar formation, 
and the observed cluster rates cannot rule this out. 
Subject headings: later 



1. INTRODUCTION 

Type la supernovae (SNe la) play a central role in 
astrophysics because of their use as standard candles 
for cosmography (e.g. lRiess et al.|[T998HPerlmutter et al.l 
[l999), as well as being a prime source for the chemical 
enrichment of the Universe. While it is widely accepted 
that SNe la are the thermonuclear explosions of carbon 
oxygen white dwarfs near the Chandrasekhar mass (e.g. 
iHillebrandt fc Niemeveii[2000f) , it is still unclear by what 
mechanism the white dwarf accretes the necessary mass, 
or if there are multiple pathways to a SN la explosion (for 
a recent review, see .Mao z ,_20_LQ,) . This knowledge about 
the SN la progenitor, and by extension the types of stel- 
lar populations that house them, will aid next generation 
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cosmological measurements since any systematic change 
in SN la properties as a function of stellar population 
age, mass or star formation history will be reflected in 
the changing galaxy population surveyed at higher red- 
shift. An understanding of the progenitor(s) of SN la 
will similarly inform the location and epoch of iron en- 
richment. 

There are currently three main progenitor/explosion 
scenarios being actively studied. The first is the sin- 
gle degenerate (SD) scenario, where a single white 
dwarf accretes material from a companion, overfilling 
its Roche lobe until it reaches the C handrasekhar mass 
(jWhelan fc Iben|[T973l:INomotoll982[ ). The double degen- 
erate (DD) scenario postulates the merger of two carbon 
oxygen white dwarfs, with a total mass which exceeds 
the Chandrasekhar limit, throug h gravitational inspiral 
and subsequen t carbon ignition (|Iben fc Tutukovl|l984l : 
iWebbinlj 119841) . Finally there is the double detonation 
scenario, a specific variant of the SD scenario, in which 
an outer He-burning shell, be ing fed from a comp anion, 
undergoes several flashes fe.g. lBildsten et al]|2007D . ulti- 
mately leading to a flnal event which ignites the underly- 
ing carbon oxygen w hite dwarf before the Chandrasekhar 
mass is reached (e.g IWooslev fc Weaver! |1994f ) - a sub- 
Chandrasekhar mass event. 

One route to distinguishing between the various SN 
la progenitor scenarios is via the delay time distribution 
(DTD). The DTD is the expected SN la rate as a func- 
tion of time following a burst of star formation, and each 
potential progenitor will have a different DTD. For in- 
stance, in the SD scenario, the DTD is driven by the 
main-sequence lifetime of the secondary star, while in 
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the DD scenario the delay time is associated with the 
time necessary for gravitational w ave emission to cause 
the orbit to decay ()Greggioll2005D . 

One approach to recovering the DTD has been through 
comparing the SN la rate in stellar populations of dif- 
ferent characteristic ages. By teasing out the relative 
SN rates in galaxies ranging from actively star forming 
to long dead, several pioneering studies found evidence 
for both 'prompt' SNe la which explode ^500 Myr af- 
ter star formation, and 'delayed' SNe la which explode 
~ 1 — 10 Gyr after explosion (TScannapieco fc B ildsten 
\2mlt. IMannucci et all [20051 [2001 : ISullivan et al.ll2006f lT 
This initially suggested that there may be two separate 
populations of SNe la (the so-called 'A-l-B' model), but 
most theoretical DTD functions predict SN la occurring 
over a broad range of a ges in any case with out invok- 
ing multiple populations ( Greggio et al.l l2008'). Much re- 
cent observational work has gone into der iving the SN 
la DTD through a variety of m ethods fe.g. iTotani et alj 
2008t IBrandt et alJ [2010: ?: IMaoz fc BadenesI [20101 



Maoz et alJ 120111: iGraur et al.l |2011[ ). Theoretically, 
both analytic (e.g. lGreggiol l2005') and binary population 
synthesis model result s (e.g. Yungc lson & Livio 2000; 
iMennekens eTall [20Tol : iRuiter et al.l l2009i [MTlt ) have 
been used to predict the functional form of the DTD 
for different progenitor models. 

A relatively clean way of probing the DTD, at least for 
stellar populations > 1 Gyr old, is by measuring the clus- 
ter SN la rate from low to high (z ~ 1.5) redshifts. Clus- 
ter red sequence (or, more loosely, early type) galaxies 
have colors and spectra which are consistent with hav- 
ing formed at high redshift (z ~ 3) in a single burst 
of star formation, followed by passive ev olution (e.g. 
[Stanford etlHUQQalETsenhardt et al.ll2008f ). Given this, 
the observed SN la rate in cluster red sequence galaxies 
as a function of time after this initial burst of star for- 
mation is a direct measure of the shape of the DTD, 
without any need to model the star formation history of 
a complex galaxy population. 

The cluster SN la rate has been measured from < z < 
1.5, although the rate uncertainties continue to be dom- 
inated by small number statistics. At z < 0.2, several 
groups measure consistent cluster rates, but have noticed 
the SN rate among cluster ellipticals is a fact or of ^2-3 
higher than in the field at the ~ 1— 2g level ( Sh aron et al.l 
[20071 : [Ma^ucci et al.l[200l iDildav e"t al.ii201(? ). If true, 
this rate enhancement may be due to a sprinkling of re- 
cent star formation in cluster ellipticals, since the SN la 
rate increas es dramatically in t he presence of star for- 
mation (e.g. ISullivan et aLll2006[ ). Cluster rate measure- 
ments at moderate to high redshift are sparser, and are 
based on < 10 SN la, not always with spectroscopic 
confirmation. These pioneering measurements include 
Hubble Space Tele scope (HST) re ii nagin g of clusters at 
0.5 < z < 0.9 by [Gal- Yam etall (|2Q0l : [Sharon et all 
(|2010| ). and results taken from a subset of the Super - 
nova Legacy Survey at z ~ 0.5 (|Graham et al.|[2008l ). 
The HST Cluster Supernova Survey has calculat ed the 
SN l a rate in clusters at 0.9 < z < 1.45 (Barb arv et al.l 
I2010D . Putting all of the above results together lead to 
constraints on the late time DTD in SN la, under the 
assumption that the bulk of cluster stars formed at high 
redshift and have evolved passively since ([Barbarv et al.l 
[2010; ,Maoz et al., i2QlQi) . Note that these earlier works 



used SN la rates over the entire cluster galaxy popula- 
tion in their DTD measurements, rather than focusing 
in on the red sequence/early type population. 

In our first paper on the Multi-Epoch Nearby Clus- 
ter Survey (MENeaCS) we described our intracluster 
SN la, and used their relative numbers to constrain the 
intracluster stella r mass fraction in our cluster sample 
(|Sand et al.ll2"oH . In this paper, we describe a new mea- 
surement of the cluster SN la rate based on 23 SNe la 
found in 57 galaxy clusters at 0.05 < z < 0.15 during the 
course of our ^2 year survey. From these measurements 
and others in the literature we constrain the late-time 
(> 2 Gyr) SN la DTD, specifically utilizing those mea- 
surements of the SN la rate in red sequence or early type 
galaxies. 

A plan of the paper follows. In §[2] we discuss the ME- 
NeaCS survey strategy and cluster sample. In § [3] we 
describe our real time image reduction pipeline for dis- 
covering and characterizing transient candidates, along 
with our final products derived from data taken over the 
entire survey. Section 4 details our spectroscopic obser- 
vations to confirm cluster SNe la. In § [5] we calculate 
the luminosity and stellar masses of the MENeaCS clus- 
ter sample and § [6] details our host galaxy photometry. 
Next, §[7] details our SN detection efficiency calculations. 
Finally, in § [S] we put all of the elements together and 
calculate the MENeaCS cluster SN la rate, along with 
our SN la rate associated with red sequence galaxies. In 
§ 9 we compare our results with both cluster and field 
SN la rates in the literature. We calculate the late time 
DTD from our own cluster red sequence SN la rate, and 
similar results in the literature, in § [TU] in order to con- 
strain the SN la progenitor. We conclude and discuss 
future directions in § [TTj Unless otherwise noted, we as- 
sume a flat ACDM cosmology with n„i—0.3, Ha— 0-7 and 
a Hubble constant of Hq = 70 km s^^ Mpc^^. 

2. MENEACS SURVEY DESIGN & GOALS 

The MENeaCS survey was designed with two primary 
scientific goals in mind: measuring the SN la rate in 
z 0.1 galaxy clusters and utilizing galaxy-galaxy lens- 
ing to measure the dark matter radius of early-type 
galaxies as a function of clustercentric distance. We de- 
tail our cluster sample and observing strategy here, with 
an emphasis on those aspects of the program most rele- 
vant for our SN science goals. 

2.1. The Cluster Sample 

To accomplish both the SN and galaxy-galaxy lensing 
goals of MENeaCS, we chose the brightest X-ray clusters 
accessible from the CFHT at 0.05 < z < 0.15. The logic 
behind this choice is simple: the most X-ray luminous 
clusters tend to have the largest stellar mass, and will 
thus have the most SNe and provide a large sample of 
galaxy lenses. The cluster sample is detailed in Table [Tj 
and consists of 57 clusters. Redshifts are those reported 
in the NASA/IPAC Extragalactic Database (NED), and 
the X-ray luminosities are those reported directly from 
the BAX online cluster database, which attempts to ho- 
mogenize luminosities in the literature and uses a Hq=50 
km s~^Mpc~^ and r2m=l cosmology. We refer the reader 
to the BAX web site for detail£3. 

^'^ http://bax.ast.obs-mip.fr/ 
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Since it will be useful to know the approximate mass 
and radius of our cluster sample when calculating SN 
rates and trends, we also list M200 and i?2oo in Table [1] 
us ing the bisector regress i on Lv — M200 scaling relation 
of IReiprich fc B ohringe? ('20021. While there are more 
up to date measurements of the slope and nor maliza- 
tion of the Lx — M relation (e.g. Stanck ct al. 120061; 
Vikhhnin et~aLl 120091: [Mantz et al.[[2"010; .Hoekstra et al l 



201 If ), our chosen relation is consistent with most lit- 
erature values at the ~lfT level, and we seek consis- 
tency with our pr evious results on the IC SNe in the 
MENeaCS survey (jSand et al.ll201lD . Ah but six of the 
MENeaCS clusters can be probed out to i?200 given the 
CFHT/Megacam field of view, and even for these clus- 
ters, we were able to cover out to ~ 0.9i?2oo- Note that a 
factor of two systematic uncertainty in our cluster masses 
will result in a '^25% systematic in our value of i?2oo- 
We discuss the effect that such systematic uncertainties 
in our derived i?200 value have on our SN la rate mea- 
surements in § [8] 

2.2. Observing Strategy 

Our goal was to obtain g' and r' images (with 2 x 120s 
per band) of each of our clusters once every month during 
the two- week queue run of CFHT/Megacam, typically 
with an image quality of ^I'.'O, in order to simultaneously 
fulfill our weak lensing requirements. These individual 
epoch images were used for SN discovery (see §1321), while 
the ultimate deep stacks over all epochs are also being 
used for programs on cluster physics and gravitational 
lensing. This strategy allows us to discover and follow SN 
la at the redshifts of our clusters for several months (see 
e.g. the template light curves in Figure [1] spanning our 
redshift range), and it also allows multiple opportunities 
for spectroscopic confirmation. 

We had a dedicated program for SN foUowup and 
confirmation. When possible, follow-up photometry of 
our cluster fields was obtained with the QOPrime imager 
(|Wilhams et al.ll2004 ) on the Steward Observatory 2.3m 
Bok telescope, for which we were generally allocated sev- 
eral nights per month. These observations were primarily 
for rejecting slow-moving objects identified as possible in- 
tracluster SN candidates, and for filling in SN light curves 
to better inform our spectroscopic pr iorities, using tech - 
niques similar to those described in iSand et al.l (|2008l ). 
Since these observations were used solely for these pur- 
poses, they are not discussed further in this work. Spec- 
troscopic follow up of SN candidates was primarily done 
at the MMT, either with the Blue Channel Spectrograph 
or Hectospec (which afforded a limited target of opportu- 
nity mode, since Hectospec is queue scheduled) . We typi- 
cally were allocated ~1 night per month for spectroscopic 
follow up during MENeaCS. Additionally, target of op- 
portunity observations with Gemini North were used for 
intracluster super nova candidates; those are presented in 
ISand et al.l muh . 

3. IMAGING REDUGTION & ANALYSIS 

Our primary data reduction goals were geared toward 
the nearly real time detection and calibration of tran- 
sients, to be discussed in § 3.1-3.3 below, the best of 
which could then be followed up spectroscopically. Ad- 
ditionally, we maintained a running deep stack image 



during the course of the survey of each cluster field specif- 
ically for identifying host galaxies and our hostless, in- 
tracluster events (§ 3.4). Ultimately, we made a set of 
final deep stack images for each cluster field based on the 
calibrated data from the CFHT (§ 3.5) for our final host 
galaxy photometry. 

3.1. Real Time Imaging Data Reduction 

As mentioned, our primary SN discovery images 
came from the squa re-degree CFHT/Megacam camera 
(|Boulade et al.l[200l . Me gacam has one ~ 2 week queue 
run per month, and we aimed for one epoch every queue 
run. Very occasionally, two epochs were taken in one 
queue run to make up for epochs missed due to weather, 
i.e. when the target field was setting. Each epoch con- 
sists of 2 X 120s images in g' and r', dithered by ^ 5' to 
fill in the chip gaps while ensuring the brightest cluster 
galaxy (BCG) is centrally positioned on one of Mega- 
cam's central chips. During the first epoch of observa- 
tion of a cluster, we took 4 x 120 s g' and r' images in 
order to have sufficient depth to determine the initial as- 
trometric solution. This first, slightly deeper epoch was 
used as the reference template for SN detection through- 
out the survey. No attempt was made, after the survey 
was complete, to use an optimal template or final deep 
stack image to recover SNe not found during the course 
of the survey. 

Bias subtracted and fiat fielded CFHT Megacam im- 
ages were provided by the CFHT-developed Eli xir data 
reduction system (jMagnier fc CuiUandr3 [200I within 
~24 hours of acquisition. The final Elixir-processed data 
for a given queue run used master flat fields, but these 
are only available after the queue run - too slow for our 
needs. We thus use the real-time data for SN detec- 
tion, preliminary photometry, and to build interim "deep 
stacks" during the survey, and the Elixir-processed data 
for our final deepstacks and photometric calibrations. 

For the real-time data, we generated obje ct catalogs 
with SExTRACTOR (jBertin fc ArnoutslllQQGD and deter- 
mine d the astrometric solution with SCAMP (jBertinl 
I2006D . matching objects with the 2MASS catalog. The 
solution was constantly refined, as every new epoch used 
all previous images. We combined each epoch's data with 
SWariS, using the lanczosS interpolation function for 
image resampling. 

3.2. Image Subtraction and Transient Detection 

As described above, each epoch results in one g' and 
one r' image, and our slightly deeper first epoch im- 
ages are used as the reference template. Briefly, the 
image with a sharper point spread function (PSF) was 
aligned, PSF-matched and flux-scaled to that of the 
poorer quality image, and the reference is subtracted 
from the new epoch. Sources on the difference image are 
detected using the IRAF task daofind. These sources 
are automatically culled to reject those which have de- 
clined/disappeared since the reference; sources which 
are not round; sources in close proximity to a negative 
source, indicating a moving object; and sources in groups 
or in rows/columns with multiple sources, indicating as- 
sociation with e.g. diffraction spikes or bad pixels. To 

version 2 . 1 5 . 7 ; [htt p : //terapix. iap . fr/soft /swarp | 
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make further automatic culls, we then measure source 
fluxes in four apertures of radius 2, 4, 6, and 8 pix- 
els, an integrated flux from the best fit PSF, and its 
^ representing goodness-of-fit. Sources are rejected if 
they are not detected simultaneously in both g' and r'; 
have g' > 23.5; do not have a constantly increasing flux 
in successively larger apertures; or have a > 50. 

After these automatic culls, triplet images are created 
for the remaining transient candidates, showing the tar- 
get image, template image and the difference. Visual re- 
view of the g' and r' triplets, and rejection only of those 
objects which are clearly spurious artifacts, leads to the 
final candidate list. These transient candidates were then 
roughly prioritized according to their likelihood of being 
cluster SN la to facilitate spectroscopic followup. We 
discuss this process in § 14.11 

3.3. Real Time Photometric Calibration 

Calibrated photometry for SN candidates in real time 
was essential information for our photometric and spec- 
troscopic follow-up campaigns. Queue run zeropoints 
were derived using all cluster fields covered by SDSS 
which were imaged in photometric conditions during that 
queue run. For real-time analysis, photometric epochs 
were simply identified by visual inspection of the CFHT 
SkyProbe plots of attenuation over time for the night. 
For epochs taken during non-photometric conditions, a 
correction for attenuation was derived by comparing star 
magnitudes to those from a frame taken under photomet- 
ric conditions. Such real-time photometric analysis was 
only moderately precise, good to within ~0. 1-0.2 mag- 
nitudes, but this was sufficient for comparison with our 
imaging follow-up at the Bok telescope, and vetting and 
exposure time estimation for spectroscopic follow-up. 

3.4. Real Time Host Identification 

As mentioned in § 13. 1[ we maintained interim deep- 
stack images of each cluster and source catalogs to facil- 
itate quick identification of SN host galaxies and likely 
cluster members. An interim host galaxy was assigned 
to each SN candida te using the method of the SNLS 
(jSullivan et al.l[2006[ ). which uses a host-SN elliptical sep- 
aration normalized by the candidate host galaxy size to 
assign hosts - the dimensionless J? par ameter described 
in Equation (1) of iSullivan et al.l ()2006[) . For details, we 
refer the re ader to the origin al work of ISullivan et al.l 
(|2006f), and iSand et al.l ()2011| ). who corrected a typo- 
graphical error in the original version. Four hostless, 
intracluster SNe, defined as those cluster SN la with no 
host within R=5, were also discovered and presented in 
ISand et al.l (gOUfc. 

3.5. Final deep stack images 

Beyond the running deep stack images constructed 
during the survey for on-the-fly host identification, we 
have also made a uniform set of deep stack images from 
our entire CFHT/Megacam data set. These stacks will 
be used for host photometry and colors, identification of 
the red sequence in each of our clusters, and for overall 
measurement of our cluster luminosities. Each of these 
measurements are key for determining and understand- 
ing our SN la rate. 

For these final deep stacks wc utilize "Elixir" processed 
Megacam images (|Magnier fc Cuillandr&.2004 ) from the 



Canadian Astronomy Data CenteiP^. The seeing of each 
input image was determined automatically with an initial 
run of the automated source extraction package SEx- 
TRACTOR. Any images with > I'.'O seeing were not in- 
cluded, although this amounted to only ~10% of the 
total. The images were resampled to a common grid 
with SCAMP and median-combined with SWarp, as 
discussed in § 3.1. These final deep stack images were 
calibrated directly to the SDSS filter system, including 
both a zeropoint and color term. Cluster fields with no 
SDSS overlap were calibrated from photometric nights 
on which an SDSS cluster was observed. 

For those fields that hosted cluster SNe, we have also 
made similar deep stack images including only epochs be- 
fore supernova discovery, specifically to measure uncon- 
taminated host colors. These SN-free stacks were gen- 
erated and calibrated as described above, with our SN 
host photometry presented in § [6l 

4. SUPERNOVA SPECTROSCOPY 

In this section, we detail our methodology for obtaining 
SN spectroscopy, and ultimately identifying cluster SN 
la - a critical ingredient for calculating the cluster SN la 
rate. 

4.1. Spectroscopic prioritization 

For two weeks in any given month during the ME- 
NeaCS program, we would obtain queue observations 
of our cluster fields from the CFHT and quickly iden- 
tify viable cluster SN candidates. Nearly all of our 
spectroscopy was scheduled classically, with roughly one 
night allocated per month. Often this night would be al- 
located during the middle of the Megacam queue run, so 
that a potentially good cluster SN la target might have 
to wait a number of weeks before being spectroscopically 
confirmed. Weather or instrument problems could also 
bump the spectroscopic confirmation of a SN to a subse- 
quent month. 

Nonetheless, our spectroscopic campaign was relatively 
robust to these scheduling and weather difficulties due to 
the relative brightness of our SNe (Figure [Ij, and the rel- 
atively smal l number of viable cluster SN la targets. As 
discussed in ISand et al.l (|2011l ). we focused our spectro- 
scopic follow up on likely cluster SN la, e.g. a newly 
discovered SN candidate had to be consistent with a SN 
la within <30 days of explosion at the cluster redshift. 
This rough magnitude and color cut changed depending 
on the cluster redshift but ranged from —0.5 < g—r < 0.8 
and 17.5 ^ g ^ 22.5. A good spectroscopic target which 
was not observed for weather or other reasons would be 
pursued in subsequent spectroscopy runs. Preference was 
given for good cluster SN candidates with a clustercentric 
distance of i? < 1.2i?200, although this was not a hard 
rule, especially if no other viable candidates were avail- 
able. We do note that any SN candidate that appeared 
to be hostless was given the highest spectroscopic prior- 
ity regardless of color or brightness, and w as generally 
sent to Gemini as a target of opportunity (jSand et al.l 

Eonl). 

We have also utilized the NASA/IPAC Extragalactic 
Database (NED) to identify known variable sources and 
objects with redshifts in the literature. This allows us, 

^* See |http: / / cadcwww.dao.nrc.ca/ 1 
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for example, to cull known variable objects (e.g. active 
galactic nuclei and variable stars) and to identify red- 
shifts of host galaxies. For our purposes, we ignore SN 
candidates that occur in host galaxies with redshifts in- 
consistent with cluster membership. 

During a typical CFHT queue run, we would iden- 
tify ~3-4 high priority cluster SN la candidates, which 
could be easily followed up with a night of 6m-class spec- 
troscopy, with time to spare for other transient spec- 
troscopy. There was a small fraction of time during 
the MENeaCS survey ('^9%) where we had successive 
weather and instrument related problems and were not 
able to perform spectroscopy on viable cluster SN la can- 
didates (see § 4.1 of Sand et al. 2011 for a discussion on 
this). We factor this into our SN rate calculations in § 8. 

4.2. Observations 

Spectroscopic identification of our viable cluster SN la 
candidates is critical for confirming their nature. We fo- 
cus on our cluster SN la here, while our cluster core col- 
lapse SNe are presented in a companion paper (Graham 
et al. in prep). 

Spectroscopic observations were obtained with one of 
three telescope/instrument combinations: Gemini North 
and the Gemini Multi Object Spectrograph (GMOS; 
iHook et al.ll2004D. the MMT and its Blue Channel Spec- 
trograph (BCS;'Schmidt et aLngSgf), and the MMT with 
Hectospec (Fabricant ct al. 2005). For our GMOS obser- 
vations (obtained solely for our IC SNe candidates) we 
used the R400 grating centered at 680 nm, along with 
the GG455 order blocking filter; see Sand et al. (2011) 
for the details of these observations. With BCS, the 300 
line grating with a central wavelength of 580 nm was al- 
ways used. Hectospec was used with the 270 line grating, 
giving wavelength coverage between '^360-800 nm, with 
one exception. For Abell2670.9_08_l (Figurel^), the Hec- 
tospec 600 line grating was used due to scheduling con- 
straints, giving wavelength coverage between ^550-750 
nm. We discuss the spectroscopic identification of this 
object in more detail in the next subsection. Standard 
CCD processing, one dimensional spectrum extraction 
and calibration were accomplished with IRAF for our 
GMOS and BCS spectra. The Hectospec spectra were 
pipeline processed at the Harvard-Smithsonian Center 
for Astrophysics. 

4.3. Supernova Identification 

To determine supernova type from our spectra, we use 
the publi cly available Supernova Identification (SNID) 
software (jBlondin fc Tonrvl [20071) . SNID correlates an 
input spectrum with a large library of template spectra, 
including SNe of all types, variable stars, luminous blue 
variables, active galactic nuclei (AGN) and galaxies. We 
use the same definition of a ' good' correla tion as previous 
works (|Miknaitis et al.ll2007l:lFolev et al.ll2009[ ). which we 
summarize briefly here while describing the basic SNID 
procedure. First, SNID places input and template spec- 
tra on the same logarithmic wavelength scale and divides 
out a pseudo-continuum so that correlations are not de- 
pendent on flux calibration errors or reddening. Also, 
since many of our SN spectra have significant host galaxy 
contamination, the removal of a continuum shape mini- 
mizes any false correlations that this may cause. A cor- 
relation redshift is calculated based on the highest peak 



in the correlation function, and a correlation parameter, 
r, is determined and defined as the ratio of the height of 
the correlation peak to the root mean square of the anti- 
symmetric component of the correlation function aro und 
the corr elation redshift fsee IBlondin fc Tonrvl [20071 for 
details). To ensure that correlations are meaningful, they 
are weighted by the overlap in ln(A) (the 'lap' parameter) 
between input and template spectrum. We consider, as 
in previous work, a 'good' correlation to have r(lap) — 
rx lap > 5 and lap = In (Ai/Ao) > 0.40. The quality of 
our best matches can be seen in Figures [2H5| 

For each candidate SN spectrum, we attempt to de- 
termine the type, subt ype, redshift and a ge by executing 
four SNID runs, as in iFolev et all ([200^ ). who we refer 
the reader to for details. A confident SN type is reported 
only if 50% of the 'good' SNID correlations were of that 
same type, along with the best-matched template. Like- 
wise, it was often impossible to determine a reliable SN 
subtype even if its type was clear, since we formally re- 
quire that 50% of the 'good' correlations (along with the 
best match) must be of that subtype. Even so, our sub- 
type determinations should be taken with a grain of salt 
given that they are based on a single spectrum, and we 
have only sparse light curves. A summary of our cluster 
SNe la are in Table [2|- no te that our f our IC SN la have 
already been presented in lSand et al.l (|2011[ ). but we in- 
clude them here as well for completeness. In Figures [2]-[8| 
we present all of our cluster SN la spectra, twenty-three 
in all. We also overplot the best-matched SNID template, 
which are listed in Table [5] as well. 

One exception for our standard supernova typing 
scheme was made for Abell2670_8_08_l. This object was 
observed with Hectospec with the 600 line grating (with 
a wavelength range of ^550-750 nm). Due to the lim- 
ited wavelength range, SNID had difficulty in converging 
to a solution. To narrow the search range, we input an 
initial guess redshift corresponding to that of the clus- 
ter {z — 0.072), and restricted the type of supernova 
searched to SNe la. This led to excellent fits to the avail- 
able data, as can be seen in Figure [71 Given this, we 
take this object to be a SN la at z = 0.076 ± 0.0036, in 
accordance with the median and dispersion of redshifts 
resulting in good fits. No subtype determination was 
attempted. 

We also inspect each SN spectrum for host galaxy ab- 
sorption and/or emission lines to pinpoint its redshift 
precisely, and we report this redshift in Table [2] when 
available. When these features are found, they agree 
with the SNID-derived redshift in all cases. 

4.4. Cluster Membership 

MENeaCS spectroscopically confirmed ^^70 SNe, and 
we must separate foreground/background objects from 
those associated with the target clusters. We start by 
plotting the histogram of the velocity difference between 
the cluster redshift and our SNe la {vsn — Vdus) in 
Figure [9] A cutoff value for cluster membership of 
Av < 3000 km s^^ is a natural fit for SN membership in 
a survey like MENeaCS. MENeaCS clusters have a veloc- 
ity dispersion of ~600-1000 km s^^ (Sand et al. 2011), 
while SN velocities can only be measured to ~1000-2000 
km s^^ due to their broad spectral features and intrinsic 
diversity (see Table [2]). As previously noted in §[2l ME- 
NeaCS becomes incomplete at ~ iZ200 due to the field of 
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view of Megacam. Thus, while we report aU of our SN 
la with Av < 3000 km in Table [2l we focus on those 
within either 1 Mpc or i?200 to calculate our SN rate in 
§ m If we were instead to place our Av cutoff value at 
5000 km s~^, we would only add three SN la into con- 
tention, all of which are at i? > 1.2i?2oo, and so no new 
objects would be added into our rate calculations. In the 
end, we have identified eleven cluster SN la within 1 Mpc 
of the cluster center, and sixteen within i?2oo- 

We note that our cut of Aw < 3000 km s^^ and Rsn < 
i?2oo is comparable with other studies in the literature. 
The intragroup SN study of iMcGee fc Baloghl (|2010l ) 
only included SNe with Rsn < Rvir and A y < 3000 km 
s~\ w hile the SDSS cluster SN results of iDildav et al.l 
(|2010l) used i?s7v <l-4 Mpc and Av < 4500km 8"^ The 
z ~ 0.5 - 0.9 cluster SN study of iSharon et"all (pOlO ) 
used a velocity cut of Av < 2600km s~^. 

5. CLUSTER STELLAR LUMINOSITY AND MASS 

In this section, we calculate the luminosity of each ME- 
NeaCS cluster and infer stellar masses - a crucial ingre- 
dient for measuring the SN rate. Given that we do not 
have a complete redshift survey to identify all cluster 
members, we measure our luminosities in a statistical 
sense, as has been done in previous cluster SN rate ana l- 
yses (e.g. ISharon et al.1[2007l IMot IBarbarv et aT1[20lol) . 

The basic approach is straightforward. Galaxy lumi- 
nosities are measured in each cluster deep stack image 
out to various clustercentric radii. Since these fields con- 
tain both cluster galaxies and foreground/background 
contamination, several background fields are used to 
gauge the typical contamination. K-corrections are de- 
termined for each cluster's redshift to determine rest 
frame properties and the statistical background is sub- 
tracted. From there, further corrections can be made for 
photometric completeness (utilizing a standard cluster 
luminosity function) and the inclusion of the intracluster 
light component. Luminosities of specific cluster galaxy 
populations can be determined (i.e. red sequence galax- 
ies), so that their SN rate can be measured directly. 

We used the stacked images of each cluster, described 
in i) 13.51 to create galaxy catalogs. We run SExtractor 
(jBertin fc Arnoutslll996[) in dual image mode, carefully 
setting the parameters so that faint galaxies were de- 
blended from bright nearby galaxies, as is common near 
the EGG. We use SExtractor MAG_AUT0 photom- 
etry for galaxy magnitudes, and a fixed 3 arcsecond cir- 
cular aperture for galaxy colors. Stars are removed from 
the catalogs by fiagging objects that lie on the stellar lo- 
cus in half-light radius {ri/2) versus magnitude. All of 
our photometry is corrected for foreground Galactic ex- 
tinction using E(B-V) values from Schlcgcl et al. (1998). 
As in previous cluster SN work, we discard galaxies in 
both the cluster frame and the background fields which 
are brighter than the identified BGG in each cluster. 

Once the catalogs are cleaned of stellar objects, a 
color magnitude diagram is made, and a red sequence 
is fit. This process will be discussed in detail by Bild- 
fell et al. (2012, i n prep), and is a simi lar procedure to 
that presented in iPimbblet et al.l (|2002[ ). We find tight 
red sequences in all of our clusters, with typical scatter 
about the linear fit of ^^0. 03-0. 04 mag for any given clus- 
ter. These measurements are important for two reasons. 
First, we do not include galaxies redder than the cluster 



red sequence (accounting for its scatter) into our cluster 
luminosity measurements (either in the background field 
or in the cluster field), as these objects, aside from the 
possibility of an extremely redde ned cluster member, ar e 
likely background galaxies (e.g. iGladders fc Yed [20001 ) . 
Second, we want to measure the red sequence SN la rate 
(§ [5]) directly, as these galaxies have a relatively pure 
stellar population of old, ~10 Gyr-old stars and are ideal 
for constraining the SN la delay time distribution (§[10]). 

For our 'background' fields, we use the four CFHT 
Legacy Survey Deep Fields, which cover a total of four 
square degrees. We downloade d the image s from the 
MegaPipe GFHTLS web pageR dGwvnl [20081 ) . using the 
"full" version containing all usable data. Galaxy catalogs 
are created using SExtractor in a fashion identical to 
that of our clusters fields. These deep images are com- 
plete to g ~ 25.5 and r ~ 25.5 as indicated by their 
galaxy number counts, which is deeper than our cluster 
deep stack fields. We only count background galaxies 
down to our chosen faint magnitude level for each clus- 
ter. We take the average luminosity per arcmin^ over the 
four background fields as our background value, and we 
take the standard deviation in luminosity in the four dif- 
ferent fields as our uncertainty, since this is an estimate 
of cosmic variations in large scale structure. 

Since only frames with PSF FWHM < 1.0 arcsec were 
included in our deep stack images, each cluster and each 
band have slightly different depths. To gauge the com- 
pleteness of each field, we plotted histograms of galaxy 
number counts versus magnitude, and did a direct com- 
parison with our deeper GFHTLS background fields. For 
the sake of simplicity, and because all of our cluster fields 
were complete at least down to g,r=23.5 mag, we only 
included galaxies brighter than this in our cluster lumi- 
nosity calculation. 

We use K-corrections derived f rom the kcorrect soft - 
ware package, as described in iBlanton fc RoweisI (|2007D . 
which we refer the reader to for details. Briefly, KCOR- 
RECT fits a spectral energy distribution (SED) to galaxy 
fiuxes based on a linear combination of a small num- 
ber of known galaxy templates, which are derived from 
the Bruzual-Gharlqt ste llar evolution synthesis code 
(jBruzual fc Gharlotl[2003[ ) ■ This SED is then used to find 
the K-correction necessary to measure absolute g and r 
band magnitudes. Additionally, since the templates can 
be interpreted physically, the SED fit also provides an 
estimate of the stellar mass-to-light ratio. We use these 
mass-to-light ratios to convert from stellar luminosity to 
stellar mass. We investigate the limitations of using only 
two photometric bands {g and r) for our K-corrections, 
by comparing with 5-band SDSS photometry when there 
is overlap with the MENeaGS sample in § 15.11 

For each cluster, we sum the K-corrected g and 7'-band 
luminosity of all galaxies brighter than our chosen de- 
tection limit of g, r=23.5 mag, fainter then the cluster's 
BGG, and within our red-color cut to exclude objects 
redder than our measured cluster red sequence. We ap- 
ply an identical cut to the four K-corrected GFHTLS 
Deep Field catalogs, and remove this statistical back- 
ground. We then make a small correction to account 
for the total luminosity of all faint cluster galaxies below 

http://www2.cadc-ccda.hia-iha.nrc- 
cnrc.gc.ca/community/CFHTLS-SG /docs / cfhtls.html 
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o ur detection limit, mimicking the procedure described 
in I Sand et al.l ()2011l ). We adopt the same faint-end slope 
(a=— 1.03 ) and normalization as that found for the Virgo 
Cluster in lTrentham fc TuUvl (|2002l ) in our cluster lumi- 
nosity function. The derived correction factors - at most 
^1 to 2% of the clusters' overall luminosity - are in- 
significant in comparison to our overall luminosity mea- 
surement uncertainties, and so we do not discuss any 
further possible ambiguities associated with our choice of 
the faint end slope or normalization here. An additional 
factor is necessary to account for the intracluster light 
component of the clusters' luminosity, as we did discover 
four intracluster S N la during the course of MENeaCS 
(|Sand et al.l 1201 1[ ). We directly adopt the intracluster 
stellar mass fraction found for the MENeaCS sample, 
//CL=0.16to;Jg, and applied this factor to obtain our 
final cluster luminosities. 

As mentioned, Kcorrect also outputs a stellar mass 
to light (M/L) ratio corresponding to the best- fit galaxy 
SED. The Kcorrect code assumes a Chabrier initial 
mass function (IMF), while the cluster SN rate results 
in the literature have all been normalized to a 'diet' 
Salpeter IMF (|Maoz et al.l[2010,) . We therefore convert 
our M/ L to a 'diet' Sa l peter IMF by using the conver- 
sions of iBernardi et al.l (j2010[) . We use these M/L ratios 
to convert our cluster luminosities to a stellar mass. The 
average stellar M/L in our clusters (dividing the total 
stellar mass by luminosity) is M/L^ = 2.4, M/Lg = 3.5 
for the cluster luminosity as a whole, and M/Lr = 2.5, 
M/Lg = 3.8 for our red sequence galaxies. 

We list our cluster luminosities and stellar masses 
within i?2oo in Table |3l We include values for the cluster 
as a whole (including galaxies below our detection limit 
and the intracluster light) and for our designated red se- 
quence galaxies only. Similar values were calculated at 
radii of 1 Mpc, 1.25i?2oo and 0.75i?2oo so that we could 
gauge the sensitivity of our SN rate measurements to 
our i?2oo values and out to a fixed cluster radius (§ H]). 
Finally, we also list i3-band luminosities in Table [H as 
many SN rates have traditionally been reported per unit 
i?-band luminosity. To convert our g and r band lu- 
minosities to jj-band, we use the filt er transformations 
presented in iBlanton fc Roweii ()2007D . 

5.1. Comparison to SDSS luminosities 

To test our luminosity measurement procedure, we 
have also derived cluster luminosities for the 43 clusters 
in the MENeaCS sample which are in the SDSS foot- 
print. We downloaded u, g, r, i, z photometry from SDSS 
Data Release ^3, choosing the model magnitudes of ex- 
tended objects, and excluding point sources. We have 
also downloaded SDSS data covering three of the four 
CFHTLS Deep Fields (the fourth is not in the SDSS 
footprint), and use these as background fields. We adopt 
a similar procedure for measuring the clusters' luminos- 
ity as described above, with two exceptions. First, we 
assume faint detection limits of gj—22.1 mag, in line 
with the SDSS depth. Second, we derive K-corrections 
using the full 5-band photometry, rather than the single 
color provided by the MENeaCS g, r photometry. 

For a fair comparison, we have recalculated our ME- 
NeaCS cluster luminosities with faint detection limits 

http://www.sdss3.org/dr8/ 



of g,r=22.2 mag, and use only the three of the four 
CFHTLS Deep Fields in the SDSS footprint as a back- 
ground. Then, in principle, the only differences in our 
derived cluster luminosities between the MENeaCS and 
SDSS data will arise from differences in our K-corrections 
derived in 5 bands versus 2, and from differences in pho- 
tometry between the SDSS pipeline and our own SEx- 
TRACTOR based approach. 

The comparison between the SDSS and MENeaCS 
data sets were relatively good, with the final cluster lumi- 
nosities agreeing at the <10% level, with the MENeaCS 
data systematically brighter in r-band and dimmer in 
g-band. We are satisfied with this measurement agree- 
ment at the ~10% level given the differences in the two 
data sets. We consider systematic uncertainties of ±10% 
in our luminosities and explore the consequences for our 
SN rate measurement in § [8l 

6. SUPERNOVA HOST GALAXIES 

While a detailed study of the host galaxy properties 
of the MENeaCs sample is beyond the scope of the cur- 
rent work, it is of interest to measure our hosts' basic 
photometric properties to divide the sample between red 
sequence versus bluer galaxies, and derive SN rates for 
each respective population. A full study of the MENeaCS 
SN host galaxy properties, including ultraviolet, optical 
and near infrared photometry, along with optical spec- 
troscopy, will be presented in a future contribution (Gra- 
ham et al. in preparation). 

First, we create SN-free deep stack images for each 
cluster field which hosted a SN, analogous to our deep 
stack combination procedure outlined in § 13.51 including 
only images taken before SN discovery. Some experi- 
mentation has shown that SN light is still visible in our 
differenced images ~12 months after SN discovery, mak- 
ing our approach conservative and appropriate. We then 
perform photometry identically to that described in § [51 
and we present our hosts' r-band magnitudes and g — r 
color (within a 3" aperture) in Tabled We also list the 
clustercentric radius of the host galaxy, both in kpc and 
scaled by i?200- 

We utilize our well-defined cluster red sequences to 
measure a color offset between the red sequence and each 
of our host galaxies. This serves as an indicator of the rel- 
ative color of the host with respect to the cluster galaxy 
population as a whole, and allows us to pick out those 
cluster SN la which were hosted by red sequence galaxies 
- we also list these color offsets in Table ID A plot of our 
hosts' red sequence offsets as a function of host magni- 
tude can be seen in Figure [TUl with the horizontal dashed 
lines indicating the median scatter in the red sequence 
over the entire MENeaCS sample. We simply define a 
host to be consistent with the cluster red sequence if its 
la color uncertainty overlaps with the median scatter in 
the red sequence. In Figure [TT] we also plot our hosts' 
red sequence offsets as a function of clustercentric radius. 
For our SN hosts within 1 Mpc of the cluster center, 6 
out of 8 are consistent with the red sequence, while 9 out 
of 13 are red sequence hosts within R2oa- In the overall 
cluster SN la sample, 10 out of the 19 cluster SN la which 
were hosted (four others were hostless intracluster SNe) 
were consistent with the red sequence. We use the above 
SN numbers when deriving our SN rates at different radii 
and in different galaxy populations in § 8. Although the 
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most central host galaxies are all on the red sequence, the 
sample size is too small to say anything definitive about 
the radial distribution of host properties in MENeaCS. 

Eleven out of the 19 SNe that were hosted were within 
one half-light radius, as reported by SExtractor. A 
'by eye' r^l^ fit to the radial distribution of SNe from 
their host (see Figure 1 of Sand et al. 2011) is excellent. 
This suggests that there is no SN detection deficiency in 
our galaxy cores (as verified by our detection efficiency 
simulations in § 7) and that the SNe "follow the light" 
of the galaxy population, at least in the sense of their 
radial distribution. More detailed work on the structure 
and morphology of our hosts will be presented in the 
future. 

We finally note that the bulk of our SN host galaxies 
consistent with the cluster red sequence are on the blue 
side of the sequence. We will defer a thorough discussion 
as to whether this might be an indication of some residual 
star formation in these hosts for a future work, where we 
will gather much more detailed information about each. 

7. SUPERNOVA DETECTION EFFICIENCY 

In order to determine SNIa rates from the MENeaCS 
survey, we must have a measure of our detection effi- 
ciency. To do this we generated a population of artificial 
SNe with realistic intrinsic properties, added them to our 
data, and ran them through the SN detection pipeline 
described in Section [31 The fraction of planted SN re- 
covered was then used to derive our efficiency, described 
below, and included in our rates calculation, as described 
later in Sectional 

7.1. The Population of Fake SNe la 

For a robust evaluation of our detection efficiencies, we 
planted artificial SNe la in clusters and epochs with a va- 
riety of properties. Since a calculation of our detection 
efhciency on every single image would be prohibitive, we 
used a representative sample of our data, constructed 
a look-up table of detection efficiencies from these, and 
then applied it to our sample as a whole for our SN rate 
calculation (§|H]). We used 18 clusters distributed evenly 
in redshift, and for each cluster chose 2-6 epochs which 
covered a range of PSF FWHMs, sky backgrounds, and 
dates. This included 74 of the 580 total epochs in ME- 
NeaCS (13%). 

We plant many artificial SNe la per field, but avoid 
overcrowding, so the total number of artificial coordi- 
nates for each cluster are set by the number of galaxies 
in the field. We used three hosting situations for the 
artificial SNe: red-sequence (RS) cluster hosts, non-RS 
hosts, and hostless. We plant a SN in half of the RS 
galaxies and 1% of the other galaxies in the field, along 
with dozens of hostless SNe. Without redshifts we can- 
not know which non-RS galaxies are cluster members, 
but for the purposes of detection efficiencies we can treat 
them as such. All artificial SNe are given the redshift of 
the cluster. 

To generate artificial SNe coordinates which "follow 
the light" in galaxies (i.e. occur more often in more lu- 
minous galaxies, and are distributed proportionally to 
their radial luminosity profiles) we used a custom rou- 
tine which, given the reference image, the positions of 
potential hosts, and the number of SNe, returned a set 
of coordinates distributed proportional to galaxy light. 



This resulted in realistic distributions of host luminosity 
and SNe la host offsets for the sample of artificial SNe. 
Hostless artificial SNe were distributed randomly, but we 
did not allow a 'hostless' artificial SN to be within 10 ef- 
fective radii of any catalog galaxy. 

Our population of artificial SNe la were given a very 
wide distribution of apparent magnitudes, bracketing our 
realistic detection sensitivities. A total of 68687 individ- 
ual transient sources were simulated. The raw instru- 
mental fiux of a artificial SNe la was calculated from its 
given apparent magnitude and our real-time photometric 
calibrations (described in Section 13. 3p . A custom rou- 
tine was used to plant artificial SNe la as a Moffat profile 
source with a PSF FWHM of that image. These altered 
images were run through the transient detection and au- 
tomatic culling procedure described in Section [31 

For recovery of artificial SNe, the visual review of de- 
tection triplets was unnecessary because only obviously 
spurious artifacts, not potential SN, are rejected at this 
stage, especially within our magnitude range of inter- 
est. Indeed, our faint end detection limit, as described in 
the next section, is dominated by our spectroscopic con- 
straints, and not our ability to detect transients down to 
even fainter magnitudes. 

7.2. Recovery Statistics of Fake SNe la 

For the rates calculation described in Section [H we re- 
quire the detection efficiency as a function magnitude: 
77(m). This is the fraction of planted objects which were 
recovered, in each apparent magnitude bin. We found 
the primary two factors that affected the shape of 7?(m) 
are position on the MegaCam field of view and the image 
PSF FWHM, as shown in Figure [T2l Other surveys find 
that epoch zeropoint or sky background affect their effi- 
ciency, but CFHT's queue-scheduling resulted in nearly 
all our data being taken in very good conditions. In 
the few cases where our sky backgrounds or zeropoints 
are worse, so was the seeing, and thus accounting for 
efficiency changes with image PSF FWHM is the only 
necessary factor. 

In processing we divide every image into 25 sections 
(5x5), and found the detection efficiencies are lowest in 
the corners, and lower along the edges, than in the center 
3x3 sections of the field. There is a simple reason for this. 
Wc SWarp all exposures of an epoch into 20x20k images 
centered on the BCG. Of our multi-point dither pattern, 
only 2 pointings are done per epoch, which results in 
shallower, single-exposure regions along some edges and 
in some corners (the initial reference image has 4 point- 
ings for complete coverage). In these single-exposure re- 
gions, bad columns and chip gaps are not filled, and the 
PSF convolution process creates larger artifacts around 
these areas, degrading our detection efficiencies at the 
field edges. The effect of this reduced 7y(m) on our rates 
is minimal as the field edges contain only a small amount 
of the cluster stellar mass. The shape of T]{m) is also al- 
tered in bad seeing. When the PSF FWHM was large, 
bright transients had better recovery statistics. Although 
transients bright enough to benefit from this are rare, we 
use the appropriate detection efficiency for image PSF 
FWHM, as described in Section [51 For reference, 84% of 
our data had a PSF FWHM better than I'.'l, while 38% 
was better than 0'.'8. 

The detection efficiency cutoff at faint and bright mag- 
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nitudes, as seen in Figure [T^ are easily explained. First, 
we imposed our spectroscopic foUowup criteria of —0.5 < 
5 - r < 0.8 and 17.5 < 5 < 22.5 into our artificial SN 
recovery statistics, leading to the faint-end cutoff in de- 
tection efficiency seen in Figure[T2l The bright-end cutoff 
in detection efficiency is due in part to image saturation 
(thus our detection efficiency for bright objects is higher 
in poor seeing) and a small bug in our SN detection soft- 
ware algorithm. An important and automated step in the 
SN detection pipeline fit a PSF to any SN candidate in 
our difference images. A candidate with a high > 50 
was deemed an artifact and rejected as a legitimate SN 
candidate. However, real, bright SN candidates were be- 
ing rejected with high ^ values, due to their high signal 
to noise ratio and the inadequacy of the PSF model being 
applied to them. Although a variable PSF fitting cut 
would have allowed us to routinely detect SNe ^ 1 mag- 
nitude brighter, our final detection efficiencies allowed us 
to discover SNe over several CFHT epochs. The essential 
task is to identify our true detection efficiencies so that 
an accurate SN rate can be calculated. 

Many surveys find a degraded efficiency in galaxy 
cores. We do see this, but only in the very centers of 
the brightest (r < 16 mag) galaxies. Thus, it affects 
only a very small fraction of our surveyed mass, and is 
accounted for by the fact that our artificial SNe were re- 
alistically planted proportional to luminosity; we need no 
additional correction for this. As verification, we found 
no deficiency in the number of actual SNe discovered in 
the central regions of our hosts, as discussed in § |6l 

8. THE CLUSTER SN RATE 

In this Section we present our final cluster SN rate 
measurements, along with a thorough discussion of our 
sources of statistical and systematic uncertainty. We 
present our SN la rate within 1 Mpc of the cluster center, 
as well as within i?200 ■ Additionally, we specifically mea- 
sure the SN la rate in red sequence galaxies at both radii, 
as well as the cluster galaxy population as a whole. The 
red sequence SN rate will be used in §[10] for constraining 
the late-time delay time distribution of SNe la. 

8.1. The SNIa Rate Calculation 

Our calculation of the SN la rate is similar in spirit to 
that presented in Sharon et al. (2007) and Barbary et 
al. (2011). The rate of Type la supernovae, Ria, is: 



Ria — 



e 



(1) 



where Nja is the number of SNe la discovered and Cspcc 
is our spectroscopic incompleteness due to weather and 
scheduling difficulties, which we have calculated to be 
0.91 (see § SIT] and Sand et al. 2011 for details). The 
cluster luminosity or stellar mass, depending on which 
you calculate the SN rate with respect to, surveyed in 
the j*'' observation is represented by Mj. The control 
time, Atj, for the j**^ of Nim images is expressed as: 



At = 



r]{m{t))dt 



(2) 



which varies with field position and image FWHM, as 
described in Section [7] The integration limits, ti and t2, 
are set to practical boundaries defined by our limit on 
color at time of discovery, g — r < 0.8. 

Since we are able to observe a given cluster SN for sev- 
eral months, it is likely for it to be re-discovered in sub- 
sequent epochs, or for it to be discovered by our pipeline 
after a second observation if the object was missed in the 
first (for example, if it occurred near a chip gap or the ob- 
serving conditions were less favorable in the first epoch). 
To account for this, when multiple images of a cluster 
are taken within < 50 days of each other, we subtract 
the probability that it was also detected previously: 



(3) 



where m(t) is the supernova light curve (magnitude as 
a function of time) , and rj is the SN detection efficiency. 



where yyj is for the j"^ epoch, and 77j_i and 77j_2 are the 
detection e fficiencies from the two previous epochs (see 
I Sharon eT al. 2007, for a similar methodology). 

As described in Section jT] the shape of 77(m) varies 
across the field. The stellar luminosity /mass, Mj in the 
denominator of Eqn [1] also varies across the field. To 
account for the spatial variation of 77(m) we determine 
three control time values for the center, edge, and corner 
regions of the image. These are each multiplied by the 
fraction of the total Mj contained in each field area and 
summed. We have assumed that the cluster luminos- 
ity (and stellar mass ) profi le is a projected NFW profile 
(| Navarro et al.l 119961 119971 ) with concentration parame- 
ter, c = 3, as has been observed for the galaxy profile in 
galaxy clusters (Lin & Mohr 2004). 

8.2. Distributions of Intrinsic SN la Parameters 

Although SNe la are used as standard candl es, there 
is div ersity among their intrinsic light curves (iPhilhpsI 
Il993f ) which we must include to properly calculate Ria. 
We use a Monte Carlo method in which the rate is calcu- 
lated over many SNIa light curves, with sub- types and 
peak magnitude drawn randomly from observed distri- 
butions (described below). 

Previous cluster SN rates have not explicitly consid- 
ered the 91bg-like (sub-luminous) and 91T-like (over- 
luminous) SNIa subtypes in their rate calculations. Dis- 
tributions of SN la subtype depend on the galaxy sample 
observed; 91bg-like SNe la happen more often in early- 
type galaxies, and 91T- like SNe la prefer late- type galax- 
ies (jHowell et al.ll2009D . Recently, IDeLalJ (|2011D have 
shown that ~16% and ~10% of all SNe la were 91bg- and 
91T-like, respectively, while these ratios change to ^^34% 
and '^3% in early type galaxies. As discussed in § 14. 3i 
MENeaCS did not have the well-sampled light curves or 
spectroscopic sequences necessary to definitively deter- 
mine SN la subtype, and so our reported SN rates tech- 
nically include all SN la. We will present our SNe rates 
both with and without consideration for sub-type, allow- 
ing for simple comparison with values in the literature. 
These rates always agree within la. 

Thus, we use the template light curves for normal, 
91bg-like (s ub-luminous), and 9 1T-like (over-luminous) 
SNe la from iNugent et al.l (|2002D . and apply the follow- 
ing to obtain a light curve in apparent magnitudes, m(t): 
the randomly chosen peak _B-band magnitude; a stretch 
value based on the established stretch-brightness correla- 
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tion (e.g. lPhillipijilQQaiPerlmutter et al.lll997[ ): the clus- 
ter's redshift; the appropriate B - to g (and r)-band SNIa 
k-correction (|Hsiao et alj 120071) : and the inverse of the 
photometric corrections for the epoch. 

For the Monte Carlo procedure, S N la subtype is d rawn 
randomly from the relative rates of iLi et aD ()2011h and 
the peak B-band magnitude is drawn from the volume- 
limited luminosity function of the appropr i ate su btype, 
which are derived from Table 3 of ILi et alj (|20Tl[) . Note 
that the iLi et alJ (|2011[ ) SN la luminosity function does 
not include a correction for host galaxy extinction, and 
by drawing from it we implicitly include the effects of 
host extinction into our SN rate measurement. For nor- 
mal SNela we use a stretch value correlated with its 
peak magnitu de, as in Mr nnak = —19.3 -|- a{l — s) where 
a 1.5 ( Asti er et al.ll2006h . Altogether, this yields a re- 
alistic distribution of SN la light curves for input to the 
Monte Carlo rates calculation. The 91bg- and 91T-like 
SNe la have their own Nugent template light curves and 
require no additional stretching. 

The cluster rates measurements of many past sur- 
veys (e.g. Mannucci et al. 2008) used a Gaussian 
distribution of SNIa p eak magnitudes, as described in 
iCappellaro et al.l (|1997l ) - we also present SN rates using 
that luminosity function for comparison purposes, which 
did not explicitly include the 91bg/91T subtypes. Note 
that we do not include an extra host galaxy extinction 
term when using this SN luminosity function for the rate 
calculation, but find r esults consisten t with the SN la 
rate derived from the ILi et al.l (j2011| ) luminosity func- 
tion - see the next section for details. 

8.3. Results 

Putting together the elements in the previous subsec- 
tions, we proceed to calculate the cluster SN rate within 1 
Mpc and i?200, both for the red sequence population and 
for the cluster as a whole. We ran 1000 trials of our SN 
rate calculation code for each scenario aforementioned, 
varying the SN la light curve, detection efficiency, cluster 
luminosity/stellar mass and SN number, Nja, according 
to their uncertainties or assigned distributions. In this 
way, our errors for detection efficiency (see error bars 
in Figure HH), for cluster luminosity/stellar mass (dis- 
cussed in Section [5]), and for Poisson statistics on Nja 
are included in the final statistical uncertainties. Our 
reported SN rate corresponds to the median of the 1000 
Monte Carlo trials, with the quoted statistical uncertain- 
ties corresponding to the 16th and 84th percentiles (the 
68% confidence interval). 

In Table [5] we present the MENeaCS SN rate measure- 
ments for the three discussed distributions of SN la lumi- 
nosity functions: 1) the Li ct al. (2011) SN la lu minosity 
function without SN91bg/91T subtypes, 2) the iLi et ahl 
(|2011h luminosity function with SN91bg/91T subtypes 
and 3) the Gaussian SN la luminosity function of 
ICappellaro et al.l (|1997[ ). which included no SN91bg/91T 
subtypes. We note that our results for these three SNIa 
luminosity functions are consistent to within our la sta- 
tistical uncertainties. We compare our r ates with those 
in the literature in § [9l and will use our ILi et al.l (l20TTh 
SN la luminosity function results, without subtype, for 
that purpose. 

Our quoted statistical uncertainties come from several 
factors: the intrinsic distribution of light curve proper- 



ties, Poisson statistics on the SN number, the uncertainty 
on the cluster luminosity, and uncertainties in our detec- 
tion efficiencies. To determine the relative contribution 
of each factor, we also run the Monte Carlo with the 
random draw deactivated for each factor in turn, and in- 
stead use only its fiducial value. From this analysis it is 
very clear that we are still dominated by the small num- 
ber of discovered SN la, with >90% of the error budget 
originating from this. 

We are aware of two possible sources of systematic un- 
certainty. The first was identified in § 15.11 where we 
found a offset of up to ^10% in our cluster luminosities 
depending on if we used SDSS versus MENeaCS pho- 
tometry. To investigate, we rescaled our SN rates, ad- 
justing our cluster luminosities/stellar masses by ±10%. 
Since the SN rate is inversely proportional to cluster lu- 
minosity/stellar mass this systematic is straightforward 
to compute, and can be seen to be smaller than our sta- 
tistical uncertainties in all cases. Our second systematic 
is due to a possible misestimation of ii2oo in our clus- 
ter sample, based on t he £y — M200 sca ling relation of 
iReiprich fc Bohringed (|2002h - see § O If this scaling 
relation is systematically off from the true Lx — M200 
relation, than this would systematically effect our i?2oo 
values. By comparing our SN rates at 1 Mpc and i?2oo 
in Table [HI one can see that the SN rate is not a strong 
function of cluster radius. Pressing the point, we assume 
a factor of two systematic offset in our cluster masses, 
which results in a ±25% systematic in our value of R2qo- 
Recalculating our SN rates at these radii yields SN rates 
which are offset by ^10%, and we include this systematic 
uncertainty into our reported results. 

Our final systematic uncertainties in Table [5] consist 
of the two factors discussed above, and the two values 
are summed in the case of our rates at i?200- Assum- 
ing we have identified all of the appropriate sources of 
systematic uncertainty, then we can safely say that our 
SN la rate measurements are dominated by our statis- 
tical uncertainty, given that our systematics are always 
subdominant. 

9. COMPARISON TO PREVIOUS WORK 

Cluster SN rates, at several redshifts, have been re- 
ported by multiple groups - we present a summary 
of the literature results in Tabl e El We have omitte d 
the high redshift measurement of iGal-Yam et al] (|2002( ). 
given that they found only one likely SN la candidate at 
z ^ 0.8. When necessary, we correct SN rate values to 
that expected for a "diet" Salpeter IMF since this is the 
dominant IMF utilized in the cluster SN rate literature. 
We note in the Table at which clustercentric radius the 
SN la rate was taken, especially if rates at multiple radii 
were reported. Since most cluster SN rate measurements 
are not tied to i?200 7 but a physical radius, we report 
our own rates within 1 Mpc in th is Table. Likewise, we 
present our rate derived from the ILi et al.l ((20TI SN la 
luminosity function, with no 91bg/91T subtypes. 

A plot of the extant cluster SN la measurements as 
a function of redshift can be seen in the left panel of 
Figure [H The MENeaCS cluster rate, plotted as the 
star, is consistent with the literature rates at < 2; < 0.2 
to within Icr, although they are systematically higher. 

In Table [6] and the right panel of Figure [T3l we present 
SN rates in cluster red sequence/early type galaxies, 
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when those values have been calculated. Again, the ME- 
NeaCS red sequence rate is consistent with other low and 
moderate redshift measurements. In fact, the only evi- 
dence that the red se quence /early type SN rate evolves 
with redshift is via the Barbarv etaLl (poToh z ~1.1 mea- 
surement, at the ~ 2cr level. A more precise measure- 
ment of the z 0.5 rate would confirm the evolution, 
and aid in constraining the DTD (see next Section). 

If we assume that red sequence galaxies are synony- 
mous with early type galaxies, there is no evidence 
for a cluster rate enhancement in the early-type SN la 
rate, as has been reported in the literature previously 
(iSharon et al.l 120071: [Mannucci et ai1l2008t iDildav et all 
120101 but see. Cooper et al. 2009), albeit at low sta- 
tistical significance. The early-type galaxy field rate, as 
reported by the SDSS-II Supernova Survey (jPildav et al.l 
l2^l20T0l ) at z ~ 0.1, was O.OTSto^al SNur - statistically 
identical to our own r-band SN rate among red sequence 
galaxies. Likewise, the 'average' SN la rate in ellipti- 
cal galaxies in t he LOSS surve y was O.OSltJJ;!]}? SNuM 
(see Table 6 of iLi et al.l 120 lit) , in agreement with our 
own red sequence rates, although this rate may include 
both field and cluster elliptical galaxies. Mannucci et 
al. 2008 did disentangle field and cluster elliptical rates 
from the iCappellaro et al.l (|1999D rate study, and found 
0.019l°;°J| SNuM in the field, which is <la lower than 
our own measurement, which we illustrate in the right 
panel of Figure |T3| Given that our own red sequence 
rates agree with field m easurements from SDSS-II and 
ICappellaro et al.l (|1999| ) , along with the combined LOSS 
rate, to the levels of our measurement uncertainties, we 
must await higher precision field and cluster measure- 
ments before any differences will be detectable. 

A simple model for the SN la rate is the 'A-l-B' repre- 
sentation, with Ria = A X M + B X M, where Al is the 
stellar m ass surveyed and M is the current star forma- 
tion r ate (|Scannapieco fc Bildstenll2005t iMannucci et al.l 
I200I . Both A and B are constants determined from 
the data. While this model is too simplistic if the true 
SN la DTD has a broad distribution, it is still useful for 
comparison with coefficient values in the literature. Re- 
cent updates to the 'A-l-B' model have suggested that 
the rate co mponent associate d with the A value may go 
like M°-^'^ (|Smith et al.l 120 lit ), although we do not con- 
sider this further for our simple comparison. If our red 
sequence galaxy sample is composed solely of old stars, 
with no recent star formation, than their SN la rate is 
a clean measurement of the 'A' coefficient. In that case, 
we find an A value of 4.1 ± 1.7 x W'^'^ SNe yr^^ Mq\ 
For comparison, the SDSS-II SN survey recently found 
A = 2.75tj]:|^ X lO"^'^ SNe yr'^ M'^ (Smith et al. 2011), 
while the SNLS foun d 5.3 ± 1.1 x IQ-^-^ SNe yr"! 
(|Sullivan et al.ll2006h - the MENeaCS value is in agree- 
ment with both of these measurements. 

We also checked for signs of a SN rate change as a 
function of cluster mass within the MENeaCS survey it- 
self by splitting our clusters into a 'high mass' and a 
'low mass' sample, dividing them at the median ME- 
NeaCS cluster mass of 6.61 X IO^^^Mq. When we cal- 
culate our red sequence SN la rate f o r each of these 
samples separately, using the ILi et al.l (|201lD luminos- 
ity function with no 91bg/91T subtypes and a limit- 



ing clustercentric radius of 1 Mpc, we get a SN rate of 
0.0372;° °249^g^j^^) to:oo34(sys) SNuM for our high mass 
sample and 0.0450tg;^^g!^(stat) tg:oo4?(sys) SNuM for our 
low mass sample. There is no discernible change in red 
sequence SN la rate as a function of cluster mass in the 
MENeaCS sample, although we are hindered by the small 
number of SNe in the high and low mass cluster bins. 

Finally, we note that even though we do not detect a 
cluster SN rate enhancement in red sequence/early type 
galaxies versus that seen in the field, it is still possible 
that some or all of these SNe belonged to a young stellar 
population, as has been suggested to be the cause of the 
cluster elliptical rate enhancement. Even galaxies that 
lie on the optical red se quence can have some recent star 
forma tion activity (e.g. iKavirai et al.ll2007l : iTrager et al.l 
I2OO8I among others), and the fact that nearly all of our 
red sequence hosts lie on the blue side of the red sequence 
indicates that more detailed study is warranted. How- 
ever, as we constrain the late-time DTD of SN la in the 
next section, we make the simplest assumption consistent 
with our current data - red sequence/early type galaxies 
are made of stars formed at high redshift {z ~ 3) in a 
single burst, and have evolved passively until the present 
epoch. 

10. CONSTRAINTS ON THE DELAY TIME DISTRIBUTION 

In this section we combine the MENeaCS cluster SN 
la rate measurement at z 0.1 with others from the lit- 
erature at < z < 1.4 to constrain the late time DTD of 
SN la. In contrast to previous work, we focus on cluster 
rates in red sequence (or early type) galaxies which we 
assume are composed of stars that formed at high red- 
shift in a single burst, and have evolved passively to this 
day. Utilizing the cluster SN rates over all galaxy types 
risks including SN la from younger stellar populations, 
weakening the power of deriving the late-time DTD from 
cluster SN measurements. Indeed, MENeaCS has spec- 
troscopically confirmed seven cluster core collapse SNe 
(Graham et al., in preparation), even though they were 
not the primary target of the survey, suggesting that star 
formation continues in our cluster fields. 

Under the assumption that all cluster red se- 
quence/early type ga laxy stars formed in a single burst 
at high redshift (e.g. iStanford et"aLlll998l among many 
others), the recovery of the DTD is straightforward - it 
is the SN rate renormalized to the stellar mass at the 
formation epoch rather than at the epoch of measure- 
ment. To account for mass lo st during stellar evo l ution , 
we take the tabulated values of lBruzual fc Charloti (|2003D 
for a Salpeter IMF, and convert to a "diet" Salpeter 
by dividi ng the tabulated mass lost by a factor of 0.7 
(see e.g. iMaoz et al.ll2010l ). Thus, the DTD is *(t) = 
m{t)RsNia(t), where m(t) is the remaining fraction of 
stellar mass at time t after the star formation burst, 
which for our chosen IMF is m(t) = 1 — mioss{t)/0-7. 

The DTDs derived from the red sequence/early type 
cluster SN la rates in the literature are plotted in Fig- 
ure 1141 A simple, and theoretically motivated, way to 
parameterize the late-time DTD is to adopt a power law, 
^'(t) oc , where t is the time since the stellar system's 
formation. We choose to leave the normalization as a 
free parameter in our fits since we have no leverage on 
the DTD at early times using the cluster SN technique. 
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For our assumed star formation epoch of = 3, we find 
a best fit value for the late time DTD scaling exponent 
of s = -1.62 ± 54, using the MPFIT package in IDL 
iMarkwardtl (|2009l ). If we change the formation epoch to 
Zf = 2.5, we find s ~ —1.50 ± 0.49. Likewise, a for- 
mation epoch oi Zf — 4.0 yields a scaling exponent of 
s = -1.78 ±0.60. 

These values are largely consistent with similar late- 
time DTD measurements in the literature, although most 
find values closer to s ~ —1. Using SN la in ellipti- 
cal galaxies from the Subaru/XMM-Newton Deep Sur- 
vev. iTotani eTall (pOOl found s = -1.08 ± 0.15, al- 
though their SNe were not spectroscopically confirmed. 
The most recent derivation of the DTD from clus t er SN 
rate measurements have come from iMaoz et al.l (|20100 
and lBarbarv et all (|2010( 1. finding power law scaling ex- 
ponents of s = —1.2 ± 0.3 and s = — 1.3llQ'4g, respec- 
tively, for an assumed star formation burst at z/ = 3. 
Both of these analyses used the cluster SN rate over all 
galaxy types, with possible contamination from SNe la 
from young stellar populations. 



10.1. Comparison with theoretical predictions 

As touched on in § 1, there are currently three ma- 
jor progenitor /explosion paths to SNe la; the SD sce- 
nario, the DD scenario and the double detonation sce- 
nario. We will compare DTD predictions for each sce- 
nario with our measurement of 4* oc where s — 
— 1.62 ± 0.54 to constrain the possible progenitor. At 
present, DTD predictions involving binary population 
synthesis have trouble matching the observed SN la 
rates (e.g. Yungclson fc Livid 120001: iRuiter et all 120091: 
iMennekens et al.ll201Q) . and so we will focus on the shape 
of the DTD rather than on its normalization, although 
this caveat should be kept in mind. 

A generic prediction of the DD scenario is a late time 
DT D with a power law form which goes like 'i> (x t"^ 
(e.g. IGreggid[2005t IRuiter et al.llMll among others), al- 
though the specific exponential slope is dependent on the 
initial binary star separation distribution. Our power law 
scaling exponent of s = — 1.62 ± 0.54 can be considered 
consistent with the DD scenario given this slight ambi- 
guity in theoretical predictions. 

In general, DTD predictions for the SD scenario pre- 
dict late-time power law slopes which are steeper than 
for the DD scenario (Greggio 2005,). In the binary pop- 
ulation synthesis calculations of lRuiter et al.l (|2009l ). the 
SD DTD looks more like a step function with the SN 
la rate quickly dropping by roughly an order magni- 
tude ^6-7 Gyr a fter star formation. The SD DTD of 
IMennekens et al] pOlO) has a similar shape, with the 
drop off age perhaps occurring a couple of Gyr earlier. 
This is not inconsistent with our observed DTD, espe- 
cially given the current uncertainty in the cluster early 
type SN rate at z ~ 0.5, and an improved measurement 
here would be a promising pat h forward. 

Recently, IRuiter et all (|2011D determined the DTD for 
the double detonation scenario, and found a power law 
falloff with 4' oc t~^. This is in general agreement with 
our own late time DTD measurement to within the un- 
certainties^ 

Finally, iPritchet etall ()2008[ ) have proposed a SN la 
DTD which is proportional to the formation rate of white 



dwarfs, and should have a power law form which goes 
like vl/ oc While this DTD is compatible with the 

trend of SN la rate versus specific star formation rate 
among the hosts in the Supernova Legacy Survey, it is 
incompatible with ou r results at the > 2a level. Indeed, 
IPritchet et al.l (|2008l ) themselves suggest that an addi- 
tional SN la channel is necessary even in the presence of 
their model. 

To conclude, our late-time DTD measurement is in 
broad agreement with recent binary population synthe- 
sis calculations for the double detonation scenario, and 
is only ~ Icr discrepant from expectations for the DD 
scenario. An improved measurement of the cluster SN 
rate at z '--^ 0.5 would help to definitively rule out DTDs 
with shapes consistent with the SD scenario. 

11. CONCLUSIONS AND FUTURE WORK 

MENeaCS surveyed 57 X-ray selected galaxy clusters 
(0.05 < z < 0.15) over a two year span and discov- 
ered 23 cluster SN la. After carefully calculating our 
SN detection efficiency, cluster luminosities and stellar 
masses, and taking into account the known dispersion in 
SN la properties, we implemented a Monte Carlo tech- 
nique for calculating our SN la rate. Our final clus- 
ter SN la rate within i?2oo (1 Mpc) is 0.042t[j °}2+0-0jo 



SNuM (0.049: 



0.014-^0.004 



SNuM), while our SN la rate 



in red sequence galaxies is 0.041 



+0.015+0.005 
-0.015-0.010 



SNuM 



(0.04lt°;;]| ^;°;°°^ S NuM), assuming the SN luminosity 



function of lLi et all ()2011h with no SN la subtypes. First 
we note that the cluster SN la rates among red sequence 
galaxies and among the cluster galaxy population as a 
whole are statistically consistent with each other, naively 
suggesting that the cluster galaxy SN la rate is not 
strongly contaminated by SN la associated with recent 
star formation. A comparison between our red sequence 
SN la rates and those reported for red/early type galaxies 
in the literature are consistent at the ~ la level, a find- 
ing predicted by iCooper et al.l ()2009|) . Although we do 
not see a cluster red sequence SN la rate enhancement 
versus the field, our results are also consistent, within 
the uncertaint ies, with the slight enhancement seen in 
previous wor kdSharon et al.l 120071: iMannucci eFall 120081: 
iDildav et al.ll2010D . We can not conclusively rule out the 
slight early type SN la rate enhancement seen in these 
earlier works. 

We gathered cluster red sequence/early type SN la 
rates from the literature, along with our own measure- 
ments, to calculate the late time (> 2 Gyr) DTD. We 
fit a power law to the extant data - \E'(t) oc - and we 
found s = —1.62 ± 0.54, assuming a brief star formation 
epoch of Zf — 3 for the red sequence stellar population, 
followed by passive evolution. This DTD is consistent 
with both the DD and double detonation scenarios for 
the SN la progenitor. An improved measurement of the 
red sequence SN la rate at z ~ 0.5 would also constrain 
the SD progenitor model, since it is predicted to have an 
order of magnitude drop-off at roughly that epoch. 

There are several future avenues of research. A com- 
panion paper to the current work will present our cluster 
core collapse SN rate, and discuss the implications for 
current star formation in the cluster environment. Addi- 
tionally, we are gathering detailed data on the host prop- 
erties of our SNe la, and of the cluster galaxies generally. 
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It is intriguing that nearly all of our red sequence hosts 
lie on the blue side of the sequence, and we are gathering 
deep optical spectroscopy to search for recent star forma- 
tion in these systems. Is it possible that the mismatch in 
normalization between theoretical late time DTDs and 
the measured SN la rates in supposedly old stellar pop- 
ulations is because of small amoun ts of star formation 
in cluster early ty pe galaxies (e.g. iKavirai et al.l 120071 : 
iTrager et al.1120081) ? 

Finally, a new cluster SN search at z ~ 0.5 of sufficient 
size to beat down the current SN rate uncertainties would 
allow for finer discernment of the late time DTD shape. 
Additionally, the discovery of intracluster SNe would pro- 
vide the first definitive measurement of the ICL fraction 
at that redshift. 
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0, 


,091 


15: 


:58: 


:20, 


,08 


+27:14:01.1 


21.24 


19.56 


2510 


Abcll 2319 


0, 


,056 


19: 


:21: 


:10, 


,20 


+43:56:43.8 


15.78 


15.61 


2350 


Abell 2409 


0, 


,148 


22: 


:00: 


:53, 


,51 


+20:58:41.8 


7.57 


9.69 


1950 


Abell 2420 


0. 


,085 


22: 


;10: 


;18, 


.60 


-12:10:12.3 


4.64 


6.67 


1760 


Abell 2426 


0. 


,098 


22: 


;14: 


;31, 


.60 


-10:22:26.6 


4.96 


7.04 


1780 


Abell 2440 


0. 


,091 


22: 


;23: 


:56, 


.94 


-01:34:59.8 


3.36 


5.33 


1630 


Abell 2443 


0. 


,108 


22: 


:26: 


:07 


.87 


+17:21:24.4 


3.22 


5.22 


1600 


Abell 2495 


0, 


,078 


22: 


:50: 


:19, 


,80 


+10:54:13.4 


2.74 


4.58 


1550 


Abcll 2597 


0, 


,085 


23: 


:25: 


:19, 


,70 


-12:07:27.7 


6.62 


8.62 


1910 


Abcll 2627 


0, 


126 


23: 


:36: 


:42, 


,06 


+23:55:29.6 


3.25 


5.29 


1600 


RXSJ2344-04 


0, 


,079 


23: 


:44: 


:18, 


,23 


-04:22:49.3 


3.62 


5.58 


1660 


Abell 2670 


0. 


,076 


23: 


:54: 


:13, 


,60 


-10:25:07.5 


2.28 


4.03 


1490 



Cluster SN rate at z ~ 0.1 
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TABLE 2 

MENeaCS cluster SN Ia 



MENEACS ID 


UT Date 


Telescope / 
Instrument 


Type 


Subtype 


z (Gal) 


z (SNID) 


Template 


Phase 


Exposure Time 

(s) 


AK*illO/i,iQ 11 no 1 
ADellz44o_i i_Uo_i 


2009- 


11- 


24.26 


MM i /riOo 


la 




0.111 


U.iiZD (U.UUD/j 


biN you 


7.5 (3.5) 


900.0 


AK^illoc^ 11 in n 
ADelloD_ii_iy_U 


2009- 


11- 


24.15 


MM i /riUo 


la 


norm 




U.Uoo4 (U.UUo/J 




ovn /i/i a\ 
Zl.iJ (14. 4j 


900.0 


AK«110/1/1Q Q nv 1 
/\ DeilZ4'-to_o_U I L 


2009- 


09- 


19.13 


MlVi L f Dy^o 


la 




nil Qc^ 
U.lioD 


n ^ inn (c\ nn7a\ 
U.iiUU l^U.UUroJ 


oiNUoeg 


9/1 n f n^ 
Z4.U (o.UJ 


yuu.U 


A K^ill 1 1 Q c; (19 1 


2009- 


07- 


23.41 


A/TA/TT /RP^ 

iVllVi ± / lJ\_/0 


la 


norm 


n iXAAA 








1 SOO 


AKoll9'^1Q ti C\A 9 
/\D011Z0iy_D_U^_Z 


2009- 


07- 


23.38 


A/TA/TT /Rr*^ 
iVllVi L j Qy^o 


la 


norm 


n r\ac\a 
U.UDUD 








Qoo n 
yuu.U 


AK^^1119Sf^ S 09 A 


2009- 


12- 


20.51 


iviivi ± j necLO 


la 


norm 




1 071 ?0 OOfll ^ 


oiN yyee 


n n fK 
u.u yo.s) 


9700 
Z 1 uu.u 


A Koii 1 f^^n Q 1 n 

/\DGll±DOU_y LO_U 


2009- 


12- 


21.45 


iviivi ± j necLO 


la 


norm 






C!Mn9or 
oi^ uzer 


no (A 7\ 


9700 

z 1 uu.u 


A Kolll Q97 1 s n 


2009- 


03- 


17.26 


iviivi ± j necLO 


la 


norm 






cKTno^ii 

olN UOQU 




ODUU.U 


Abelll991_5_23_0 


2009- 


04- 


27.23 


MMT/Heeto 


la 


norm 




0580 (() 0029^ 


SN89B 


31 n 3 3") 


3600 


Abelll991.8.13.1 


2009- 


07- 


16.15 


MMT/Hecto 


la 


norm 




0.0526 (0.0034) 


SN94D 


29.7 (11.5) 


2700.0 


Abell2409_13.14_0 


2009- 


12- 


20.06 


MMT/Heeto 


la 


norm 




0.1458 (0.0068) 


SN02er 


2.2 (6.2) 


3600.0 


Abell795.7.13.2 


2009- 


03- 


16.17 


MMT/Heeto 


la 




0.1414 


0.1388 (0.0079) 


SN03eg 


28.4 (9.7) 


3600.0 


Abell795.9.13.0 


2009- 


11- 


24.40 


MMT/BCS 


la 


norm 


0.1363 


0.1348 (0.0027) 


SN03eg 


31.1 (7.4) 


1200.0 


ZwC10628.4_13.0 


2009- 


01- 


29.14 


MMT/Heeto 


la 


norm 




0.0822 (0.0052) 


SN03eg 


-0.3 (5.5) 


3600.0 


Abell2670.9.08.1 + 


2009- 


12- 


21.06 


MMT/Heeto 


la 






0.076 (0.0036) 


SN04eo 


57.7 (21.3) 


3600.0 


Abell2033.5.19.0 


2009- 


05- 


03.22 


MMT/BCS 


la 


norm 




0.0837 (0.0039) 


SNOSdu 


20.2 (23.4) 


1200.0 


MKW3S.4.14.0 


2009- 


05- 


03.43 


MMT/BCS 


la 


norm 




0.0526 (0.0070) 


SN94ae 


16.5 (75.8) 


1200.0 


Abell21_4_19_0 


2008- 


11- 


05.07 


MMT/BCS 


la 


norm 


0.0947 


0.0913 (0.0110) 


SN99ee 


7.9 (8.7) 


1200.0 


Abell644_5.19.0 


2008- 


12- 


19.32 


MMT/BCS 


la 


norm 




0.0656 (0.0036) 


SN89B 


25.4 (17.2) 


2400.0 


RXCJ0736p39.2.09.0 


2008- 


12- 


29.28 


MMT/BCS 


la 


norm 


0.1225 


0.1236 (0.0035) 


SN94D 


18.9 (9.4) 


1200.0 


Abell85.6.08.0 t 


2009- 


07- 


04.60 


Gemini/GMOS 


la 


Ia-91bg 




0.0617 (0.0007) 


SN91bg 


34.1 (5.5) 


1800.0 


Abell2495.5.13.0 t 


2009- 


06- 


18.58 


Gemini/GMOS 


la 


norm 




0.0796 (0.0032) 


SN98aq 


83.3 (15.1) 


3000.0 


Abell399.3.14.0 t 


2008- 


11- 


28.49 


Gemini/GMOS 


la 


norm 




0.0613 (0.0025) 


SNOSdu 


14.4 (2.4) 


1200.0 



t Intracluster SN la. as presented in lSand et al.l (|2011f ) 
Values in parentheses are the standard deviation as reported by SNID. 
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Sand et al. 



TABLE 3 

Cluster Luminosities & Stellar Masses within _R2oo 



g-band Luminosity (lO^^L©) r-band Luminosity (lO^^L©) B-band Luminosity (lO^^L©) Stellar Mass (lO^^M©) 



Cluster 


Total 




RS 


Total 




RS 


Total 




RS 


Total 


xlo 


Abell2703 


1.67 ±0.51 


1, 


.00 ± 0.02 


2.54 ± 0.76 


1, 


,53 ± 0.04 


1.42 ±0.43 


0, 


,84 ±0.02 


6.19 ± 1.71 


3.94 ± 0.12 


Abcll7 


3.61 ± 0.64 


1 


.90 ± 0.04 


5.08 ± 0.94 


2, 


.76 ±0.09 


3.16 ± 0.59 


1, 


,65 ±0.04 


12.0 ± 2.07 


6.81 ± 0.24 


Abcll21 


1.89 ± 0.45 


1 


.20 ± 0.03 


2.67 ±0.60 


1, 


.78 ±0.07 


1.66 ± 0.40 


1, 


,04 ± 0.03 


6.38 ± 1.31 


4.44 ± 0.18 


Abell85 


3.16 ± 0.55 


2 


.10 ±0.02 


4.78 ± 0.95 


3, 


.09 ±0.10 


2.70 ± 0.50 


1, 


,81 ±0.03 


11.2 ± 2.05 


7.54 ± 0.26 


AbcUlig 


2.51 ± 0.29 


2 


.12 ± 0.02 


3.87 ± 0.59 


3, 


.03 ±0.07 


2.13 ± 0.27 


1, 


,85 ±0.02 


9.20 ± 1.25 


7.30 ± 0.17 


Abciuaa 


1.49 ± 0.40 


1 


.24 ±0.01 


2.35 ± 0.65 


1, 


.85 ±0.07 


1.25 ± 0.34 


1, 


,06 ±0.02 


5.67 ± 1.41 


4.57 ± 0.18 


RXCJ0132m08 


1.77 ± 0.79 





.74 ± 0.04 


2.46 ± 1.12 


1, 


.16 ±0.08 


1.56 ± 0.66 


0, 


,62 ±0.04 


5.62 ± 2.52 


3.00 ± 0.21 


Abcll399 


4.90 ± 0.48 


3 


.40 ± 0.07 


7.33 ± 0.93 


4, 


.97 ±0.11 


4.20 ± 0.47 


2, 


,94 ±0.07 


17.4 ± 2.02 


12.1 ± 0.28 


RXCJ0352pl9 


2.04 ±0.62 


0, 


.77 ±0.05 


2.67 ±0.92 


1, 


.13 ±0.11 


1.84 ±0.57 


0, 


,67 ±0.05 


5.75 ± 2.01 


2.76 ± 0.28 


Abell478 


8.57 ± 1.04 


2 


.64 ± 0.07 


10.5 ± 1.65 


3, 


.76 ± 0.20 


7.92 ± 1.07 


2, 


.31 ± 0.08 


21.7 ±3.51 


8.90 ± 0.49 


Abell553 


1.93 ±0.27 


0.83 ±0.01 


2.53 ± 0.41 


1, 


.21 ± 0.03 


1.74 ±0.26 


0.72 ± 0.01 


5.63 ± 0.87 


2.97 ± 0.08 


ZwC10628 


8.20 ± 0.49 


2.17 ±0.03 


9.43 ± 0.78 


3.08 ± 0.09 


7.79 ± 0.55 


1.90 ± 0.03 


19.1 ± 1.70 


7.39 ± 0.24 


RXCJ0736p39 


1.66 ± 0.80 


0.90 ± 0.04 


2.50 ± 1.14 


1, 


.38 ± 0.07 


1.42 ±0.63 


0.76 ± 0.03 


5.99 ± 2.53 


3.56 ± 0.20 


Abell644 


4.29 ± 0.72 


2 


.00 ± 0.03 


6.12 ± 1.17 


3, 


,02 ±0.11 


3.75 ± 0.68 


1, 


,70 ± 0.04 


14.1 ± 2.55 


7.48 ± 0.30 


Abcll646 


2.85 ± 0.75 


1 


.35 ± 0.04 


4.19 ± 1.09 


2, 


,13 ±0.07 


2.46 ± 0.65 


1, 


,13 ±0.04 


10.1 ± 2.45 


5.54 ± 0.20 


Abell655 


3.95 ± 0.72 


2 


.40 ± 0.04 


5.98 ± 1.11 


3, 


,72 ±0.07 


3.37 ± 0.64 


2, 


,02 ±0.04 


14.7 ± 2.49 


9.63 ± 0.20 


Abell754 


3.83 ± 0.32 


3 


.11 ± 0.07 


6.29 ± 0.98 


4, 


,64 ±0.11 


3.16 ± 0.33 


2, 


,67 ±0.07 


15.0 ± 2.12 


11.3 ± 0.27 


Abell763 


0.69 ± 0.47 





.45 ± 0.05 


1.06 ± 0.70 


0, 


,66 ±0.09 


0.58 ± 0.35 


0, 


,39 ±0.04 


2.42 ± 1.53 


1.61 ± 0.24 


Abcll780 


1.52 ± 0.48 





.91 ± 0.01 


2.33 ± 0.81 


1, 


,35 ±0.08 


1.29 ± 0.40 


0, 


,78 ±0.02 


5.24 ± 1.74 


3.25 ± 0.20 


Abcll795 


4.51 ± 1.01 


2 


.61 ± 0.04 


6.99 ± 1.43 


4, 


,17 ±0.09 


3.81 ± 0.87 


2, 


,18 ±0.04 


17.7 ± 3.22 


10.9 ± 0.24 


Abcll961 


2.43 ± 0.67 


1 


.75 ± 0.04 


3.84 ± 0.98 


2, 


,71 ±0.07 


2.04 ± 0.56 


1, 


,48 ±0.04 


9.71 ± 2.19 


6.99 ± 0.20 


Abell990 


3.04 ± 1.05 


2, 


.00 ± 0.05 


4.63 ± 1.47 


3, 


,12 ±0.10 


2.59 ± 0.87 


1, 


,69 ±0.05 


11.4 ± 3.30 


8.10 ± 0.27 


ZwC11023 


2.20 ± 1.04 





.99 ± 0.06 


3.25 ± 1.58 


1, 


.56 ±0.12 


1.90 ± 0.84 


0, 


.84 ± 0.06 


7.47 ± 3.57 


4.02 ± 0.33 


Abelll033 


2.88 ± 0.83 


1, 


.69 ± 0.04 


4.44 ± 1.17 


2, 


.63 ± 0.08 


2.44 ± 0.70 


1, 


.43 ± 0.04 


11.1 ± 2.60 


6.83 ±0.23 


Abelll068 


3.30 ± 1.06 


1, 


.03 ± 0.05 


4.75 ± 1.57 


1, 


.60 ± 0.09 


2.87 ±0.92 


0.87 ± 0.04 


11.1 ± 3.53 


4.14 ±0.26 


Abellll32 


3.81 ± 0.91 


2 


.36 ± 0.04 


5.78 ± 1.30 


3, 


,73 ± 0.08 


3.25 ± 0.79 


1, 


,98 ± 0.04 


14.3 ± 2.92 


9.74 ±0.21 


Abelll285 


3.62 ± 0.71 


1 


.88 ± 0.04 


5.43 ± 1.16 


2, 


,87 ±0.07 


3.10 ± 0.64 


1, 


,60 ±0.04 


13.1 ± 2.59 


7.34 ± 0.20 


Abelll348 


3.05 ± 0.73 


1 


.20 ±0.05 


4.55 ± 1.09 


1, 


,84 ±0.09 


2.62 ± 0.64 


1, 


,02 ±0.04 


11.0 ± 2.43 


4.73 ± 0.24 


Abelll361 


1.52 ± 0.78 





.83 ± 0.04 


2.19 ± 1.16 


1, 


,25 ±0.08 


1.32 ± 0.63 


0, 


,71 ±0.04 


4.93 ± 2.59 


3.18 ± 0.22 


Abelll413 


6.18 ± 1.41 


2 


.65 ± 0.08 


8.83 ± 2.08 


4, 


,09 ±0.15 


5.39 ± 1.27 


2, 


,25 ±0.08 


21.0 ± 4.67 


10.5 ± 0.41 


ZwC11215 


2.01 ± 0.66 


1 


.70 ± 0.04 


3.31 ± 1.05 


2, 


,59 ±0.10 


1.65 ± 0.54 


1, 


,44 ±0.04 


8.24 ± 2.30 


6.54 ± 0.28 


Abclll650 


2.59 ± 0.75 


1 


.91 ± 0.05 


4.21 ± 1.26 


2, 


,85 ±0.12 


2.15 ± 0.62 


1, 


,64 ±0.05 


10.1 ± 2.76 


7.05 ± 0.31 


Abclll651 


3.43 ± 0.81 


2 


.24 ± 0.04 


4.96 ± 1.25 


3, 


,33 ±0.13 


2.98 ± 0.72 


1, 


,93 ±0.04 


11.6 ± 2.71 


8.23 ± 0.34 


Abclll781 


1.69 ± 0.50 





.85 ± 0.02 


2.50 ±0.79 


1, 


,28 ±0.08 


1.45 ± 0.43 


0, 


,72 ±0.02 


5.78 ± 1.71 


3.12 ± 0.20 


Abelll795 


3.52 ± 0.71 


2, 


.02 ± 0.05 


4.98 ± 1.19 


2, 


,91 ±0.13 


3.08 ± 0.66 


1, 


,75 ±0.06 


11.0 ± 2.55 


6.98 ± 0.32 


Abelll927 


2.58 ± 0.52 


1, 


.33 ± 0.04 


3.58 ± 0.78 


2, 


.00 ± 0.08 


2.28 ± 0.48 


1, 


.14 ±0.04 


8.29 ± 1.72 


5.05 ±0.23 


Abelll991 


1.20 ±0.27 


0.80 ± 0.02 


1.83 ± 0.47 


1, 


.18 ± 0.04 


1.02 ± 0.24 


0.68 ± 0.02 


4.24 ± 1.03 


2.86 ±0.11 


Abell2029 


7.12 ± 1.20 


3 


.31 ±0.09 


10.2 ± 1.95 


4.96 ± 0.19 


6.20 ± 1.12 


2, 


.84 ± 0.09 


23.9 ± 4.29 


12.3 ±0.51 


Abell2033 


3.19 ±0.52 


1, 


.88 ± 0.04 


4.61 ±0.84 


2, 


.80 ± 0.09 


2.77 ±0.48 


1, 


,61 ±0.04 


10.8 ± 1.83 


6.97 ±0.24 


Abell2050 


2.92 ± 0.55 


1, 


.59 ±0.02 


4.28 ± 0.78 


2, 


,44 ±0.04 


2.52 ± 0.50 


1, 


,35 ± 0.02 


10.4 ± 1.76 


6.29 ± 0.13 


Abell2055 


3.72 ± 0.60 


1 


.85 ± 0.03 


5.25 ± 0.89 


2, 


,78 ±0.07 


3.26 ± 0.56 


1, 


,58 ±0.03 


12.4 ± 1.94 


6.99 ± 0.18 


Abell2064 


1.40 ± 0.50 


1 


.27 ± 0.04 


2.42 ± 0.88 


1, 


,90 ±0.07 


1.13 ± 0.39 


1, 


,09 ±0.04 


5.84 ± 1.95 


4.72 ± 0.20 


MKW3S 


0.86 ± 0.25 





.75 ± 0.01 


1.52 ± 0.51 


1, 


,13 ±0.05 


0.69 ± 0.20 


0, 


,64 ±0.02 


3.58 ± 1.11 


2.74 ± 0.13 


Abell2065 


3.66 ±0.46 


2 


.40 ± 0.03 


5.30 ±0.68 


3, 


,56 ±0.06 


3.18 ±0.43 


2, 


,07 ±0.03 


12.7 ± 1.47 


8.87 ± 0.17 


Abcll2069 


5.06 ±0.74 


3 


.10 ±0.04 


7.33 ± 1.09 


4, 


,67 ±0.08 


4.39 ± 0.68 


2, 


,65 ± 0.04 


17.8 ±2.43 


11.8 ±0.23 


Abcll2142 


7.06 ± 1.20 


3 


.98 ± 0.06 


10.3 ± 1.81 


6, 


,07 ±0.15 


6.09 ± 1.09 


3, 


,38 ±0.06 


25.2 ±4.00 


15.4 ±0.41 


Abcll2319 


15.9 ±0.62 


4 


.69 ± 0.02 


19.4 ± 1.04 


6, 


,53 ±0.09 


14.8 ±0.68 


4, 


,13 ±0.03 


41.4 ±2.24 


15.5 ± 0.22 


Abell2409 


5.13 ± 1.25 


2, 


.51 ± 0.08 


7.50 ± 1.75 


3, 


,90 ±0.14 


4.43 ± 1.10 


2, 


,12 ±0.07 


18.3 ±3.93 


10.1 ± 0.37 


Abell2420 


2.62 ±0.62 


1, 


.53 ± 0.04 


3.94 ±0.98 


2, 


,31 ±0.11 


2.24 ±0.54 


1, 


,31 ±0.05 


9.37 ±2.15 


5.76 ±0.28 


Abell2426 


4.77 ± 0.73 


2, 


.47 ± 0.04 


6.64 ±1.12 


3, 


,63 ±0.11 


4.20 ± 0.69 


2, 


,13 ±0.05 


15.4 ± 2.46 


9.04 ± 0.30 


Abell2440 


1.95 ± 0.57 


1, 


.70 ± 0.03 


3.20 ± 0.93 


2, 


.57 ± 0.08 


1.61 ±0.47 


1, 


.45 ± 0.03 


7.96 ± 2.04 


6.49 ± 0.21 


Abell2443 


2.81 ± 0.56 


1, 


.66 ± 0.02 


4.31 ± 0.89 


2, 


.53 ± 0.06 


2.39 ± 0.50 


1, 


.41 ± 0.02 


10.6 ± 2.00 


6.46 ±0.17 


Abell2495 


1.59 ± 0.47 





.90 ±0.03 


2.41 ± 0.76 


1, 


.38 ± 0.09 


1.36 ± 0.40 


0, 


.77 ± 0.03 


5.63 ± 1.67 


3.42 ±0.23 


Abell2597 


1.64 ±0.81 


1, 


.14 ±0.04 


2.72 ±1.31 


1, 


,74 ±0.10 


1.35 ±0.61 


0, 


,97 ±0.04 


6.51 ±2.91 


4.37 ±0.28 


Abell2627 


2.87 ±0.58 


1 


.17 ± 0.03 


3.97 ±0.86 


1, 


,79 ± 0.06 


2.53 ±0.54 


0, 


,99 ±0.03 


9.25 ± 1.92 


4.59 ±0.16 


RXSJ2344m04 


1.52 ±0.60 


1 


.19 ±0.04 


2.51 ±0.97 


1, 


,82 ±0.08 


1.25 ±0.48 


1, 


,02 ±0.04 


6.05 ± 2.12 


4.55 ± 0.22 


Abell2670 


3.41 ± 0.44 


2, 


.38 ± 0.02 


5.15 ±0.76 


3, 


.56 ± 0.07 


2.91 ±0.41 


2, 


.04 ± 0.02 


12.4 ± 1.65 


8.85 ±0.18 



Cluster SN rate at 2; ^ 0.1 
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TABLE 4 
MENeaCS SN Ia host properties 



MENEACS ID 


a 


S 


r 




g — r 


A(g — r)rs 




D / F?, 




(J2000.0) 


(J2000.0) 


(mag) 




(mag) 


(mag) 


(kpc) 




Abell2443.11.08.1 


22;25;59.56 


+ 17:13:32.8 


16.70 + 0, 


.02 


1.00 + 0.04 


-0.02 


960 


0.60 


Abell85_ll_19_0 


00;41;12.91 


-09:09:00.3 


15.88 + 0, 


.02 


0.87 + 0.04 


-0.03 


840 


0.53 


Abell2443.8.07.1 


22;27;07.87 


+ 17:14:36.6 


17.48 + 0, 


.02 


0.83 + 0.04 


-0.17 


1880 


1.17 


Abellll9_5.02_l 


00:56:46.57 


-01:35:47.4 


15.12 + 0, 


.02 
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TABLE 5 
MENeaCS Cluster SN Ia Rates 



Gal Rad Ns n 
Type^ SNugt' 


SN Rate 
SNur-^ SNuBd 


SNuM'^ 


With the SNIa LF of Li et al. (2010), no subtypes. 



1 Mpc b U.i5,i_,j 055_(, U.iUd.Q 037-0.009 U. iSU_o.o6.5 -0.016 U-'J4i_„.oi.s _o.004 

All 1 ^/r,^^ 11 n ifis+0'058 +0.019 r, 7+0.039+0.013 p, , qo+o. 067 +0.021 n/iQ+0-"i6+o.oo.'i 

^ '■^^P'^ "^-J^O^-0. 048 -0.015 ^-J^J^'-0. 034-0. oil '^•^^^•^-0. 054 -0.018 '^-'^^-'-0. 014 -0.004 

pq R Q n 1 i;^+0 0B8 +0.017 n ipo+0. 038+0. Oil n lon+0. 066 +0.020 r, 0/11 +0 015 +0.005 

"■^ "-200 y '-'■103-0.054-0.036 "J-l^-^-0. 036-0.023 ^'•^""-0. 063-0.041 '-'■'^^l-0. 015 -0.010 

All Ifi n 144+0 039 +0-044 f. 1 nn+0 027+0.026 p, , -,.+0.048 +0.056 049+0 012 +0.010 

^" K20O to U-144_o 035-0.034 ^"-0.024 -0.021 U-lP4_o,o4i _o,o39 U-U4i_o oiQ -O.008 

With the SNIa LF of Li et al. (2010), including subtypes. 

""^^i 1 Mr,r fi n 104+0. 1^+0.022 „ , ,_+U. 103 +0.U14 „ ^+0.166+0.024 n nt;9+0-0y« +0.006 

J. ivipo o ^-1^^-0.086-0.018 ^-1^' -0.056-0.012 "-^1^^-0.092-0.020 "•"'-'^-0.024-0.005 

All 1 Mnf 11 f, 170+0.094+0.019 p 1,4+0.083+0.014 +0.131 +0.023 „ f,p-f,+0.024 +0.006 

■"^11 1 ™P'' 11 "-l''^-0. 053-0. 016 "-1^^-0.040-0.011 "-^"^-0.063-0.019 "-""^"-0. 015 -0.005 

pee R„„„ Q r, 1SQ+0-127+0.021 r, , rjc+0. 094 +0.014 p, r,, r,+0. 151 +0.024 r, p,.Q+0. 041 +0.005 

1^'-' 11-200 » "-lo^-0. 074-0. 044 "■ l^'^-0. 048 -0.028 "-^1^-0.079-0.050 "-"^-"-0. 018 -0.010 

All R„„„ 1« r, 14Q+0-068 +0.046 r, 1 nfi+0-061 +0.030 r, 1 77+0-098 +0.056 r, 04 ■^+0-020 +0.012 

1^200 1" "-1^^-0.041-0.034 "- 1""-0. 029 -0.022 "-^ -0.051 -0.046 "-"^'-'-0.011-0.008 

With the SN la LI' ol Cappcllaro et al. (1999). 

TTI^ R n ic:o+0.0«0+U.01« n 10/1+0.053+0.012 n i oQ+0.Uy4 +0.U2U n n/1O+0-021 +0.00b 

"-^ 1 IVlPC f U.it)»_o. 058-0. 014 U.iU4_o o37_o.009 U.i«.3_o.o66 -0.017 ^-^^^-0. 015 -0.004 
All 1 Mr.f 11 f, 171+0-066+0.019 /. i r,r,+0. 041 +0.014 qq+O. 070 +0.022 „ „^„+0. 018 +0.006 
■^'1 1 11 "-I'l-0. 049-0. 016 "-1^^-0.036-0.011 "-1^^-0.059-0.018 "-"^^"-0. 015 -0.005 
Roe ■R„„„ Q p, ,f.r,+0. 056+0. 018 p ,p7+0. 040 +0.012 p, , oq+0. 067 +0.021 p. p. .0+0.015 +0.005 
1^200 9 U.lbO_(, 0.55-0.037 lU7_o o37 -0.023 U.189_o.o65 -0.044 *J-U4d_o.oi(i -O.OlO 
A]] R„„„ 1« 148+0-044+0.049 „ .04+0-028+0.027 170+0-049+0.059 o n44+°-°l^ +0-0"' 
fi-200 IP "-l^P-0. 035 -0.034 "■l"^-0. 026-0. 023 "-1 '"-0.041 -0.041 "-"^^-0.011-0.009 



Note — The first and second uncertainty intervals are the Statistical and Systematic errors, respectively. 
^ Wc present our SN la ratca in both our red sequence (RS) galaxies and Over all cluster galaxies (All). 
^ SNug = SNe(100 yr Ifll" ig,©)"! 
^ SNur = SNe(100 yr IqI'^ L^ q)~'^ 

"1 SNuB = SNe(100 yr IqI" ifl,©)"! 
^ SNuM = SNe(100 yr IqI" Mq)-1 
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TABLE 6 

Summary of Cluster SN Rate Measurements 



Reference 


Rcdshift 


SN Ratet 




(Range) 


SNuM 


Overall Cluster Rate 



This work ^ 
Mannucci et al. 2008*^ 
Dilday et al. 2010 
Sharon et al. 2007 
Dilday et al. 2010 
Graham et al. 2008 " 
Sharon et al. 2010'' 
Barbary et al. 2010 



0.10 (0.05 - 0.15) 
0.020 (0.005 - 0.04) 
0.084 (0.03 - 0.17) 
0.15 (0.06 - 0.19) 
0.225 (0.100 - 0.300) 
0.46 (0.20 - 0.60) 
0.60 (0.50 - 0.89) 
1.12 (0.90 - 1.45) 



0.057^ 



049+" "^^ 

f0.021 
-0.016 

0.0601^2? 
0.098l«:0g 

n 1 '7'7+0.212 
'-0.124 
Q it-1+0.138 
U.101_o.il6 

'^■■^"-0.19 



Cluster early-type galaxies or red sequence galaxies 



041+" "^* 
066+'' °^'' 

U.UDD_g Q20 

n5i+'' °3i 

U.UOi_Q 

0.080tHl 
0.085t[J:go? 

7+0.22 
-0.18 



This work^ 

Mannucci et al. 2008 ^ 
Dilday et al. 2010'= 
Dilday et al. 2010"^ 
Graham et al. 2008 
Barbary et al. 2010 ^ 



0.10 (0.05 - 0.15) 
0.020 (0.005 - 0.04) 
0.084 (0.03 - 0.17) 
0.225 (0.100 - 0.300) 
0.46 (0.20 - 0.60) 
1.12 (0.90 - 1.45) 



0.37^ 



Quoted SN rate 
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alu 
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Fig. 1. — Template SN ia light curves at z = 0.05 and z = 0.15 l INugent et al.l l20021. the redshift range of the MENeaCS survey. We have 
a self-imposed detection limit for SNe of g ~ 22.5 mag to accommodate efficient spectroscopic foUowup, as indicated by the dot-dashed line 
in the figure. We are capable of spectroscopically identifying our cluster SNe la roughly two or more months after explosion. 
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Fig. 2. — Spectra of cluster SN la, along with the best-fitting template SN as determined by SNID. See S I4.3l for details of the spectroscopic 
typing and Table[2]for results. SNID removes a pseudo-continuum from each spectrum before fitting to minimize host galaxy contamination. 
Top left - Abelll991.8.13.1, a cluster normal SN la at 2 = 0.0526. Top right - Abelll991.5.23.0, a cluster normal SN la at z = 0.0580. 
Bottom left - Abell2409.13.14_0 is a normal cluster SN la at z = 0.1458. Bottom right - Abelll927.5.18.0 is a normal cluster SN la at 
z = 0.0945. 
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Fig. 3. — Spectra of cluster SN la, along with the best-fitting template SN as determined by SNID. See ii l4.3l for details of the spectroscopic 
typing and Table[2]for results. SNID removes a pseudo-continuum from each spectrum before fitting to minimize host galaxy contamination. 
Top left - Abelll285.8.02_4, a cluster normal SN la at z = 0.1071. Top right - ZwC10628.4_13_0 is a normal cluster SN la at z = 0.0822. 
Bottom left - Abell795.7.13.2 is a cluster SN la at z = 0.1414 Bottom right - Abell2443.11.08.1 is a cluster SN la at z = 0.111. 
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Fig. 4. — Spectra of cluster SN la, along with the best-fitting template SN as determined by SNID. See S I4.3l for details of the spectroscopic 
typing and Table[2]for results. SNID removes a pseudo-continuum from each spectrum before fitting to minimize host galaxy contamination. 
Top left - Abell85.11.19., a normal cluster SN la at z = 0.0534. Top right - Abell2443_8.07.1 is a cluster SN la at 2 = 0.1135. Bottom 
left - Abell2319.6.04_2 is a normal cluster SN la at 2 = 0.0606. Bottom right - Abellll9.5.02_l is a normal cluster SN la at z = 0.0444. 
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Fig. 5. — Spectra of cluster SN la, along with the best-fitting template SN as determined by SNID. See S I4.3l for details of the spectroscopic 
typing and Table[2]for results. SNID removes a pseudo-continuum from each spectrum before fitting to minimize host galaxy contamination. 
Top left - Abell2033.5.19.0, a normal cluster SN la at z = 0.0837. MKW3S.4_14_0 is a cluster SN la at z = 0.0526. Abell21_4_19.0 is a 
normal cluster SN la at z = 0.0947. Abell644_5.19.0 is a normal cluster SN la at z = 0.0656. 
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Fig. 6. — A spectrum of a cluster SN la, along with the best-fitting template SN as determined by SNID. See ij 14.31 for details of the 
spectroscopic typing and Table|2]for results. SNID removes a pseudo-continuum from each spectrum before fitting to minimize host galaxy 
contamination. RXCJ0736p39.2_09_0 is a normal cluster SN la at z = 0.1225. 
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Fig. 7. — A spectrum of a cluster SN la, along with the best-fitting template SN as determined by SNID. See ij 14.31 for details of the 
spectroscopic typing and Table|2]for results. SNID removes a pseudo-continuum from each spectrum before fitting to minimize host galaxy 
contamination. Due to Hectospec scheduling, it was necessary to get this spectrum in 600 line mode, only allowing for wavelength coverage 
between ~5500 to 7500 A. Nonetheless, it is a cluster SN la, although its subtype was not obtainable due to the wavelength coverage. 
Abell2670_9.08.1 is a cluster SN la at 2 = 0.076. 
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Fig. 8. — The four intracluster SN la found by MENeaCS, and first reported in lSand et al.l l|20111 'l. See that work for details. Top left - 
Abelll650_9.13.0, a cluster normal SN la at z = 0.0.0836. Top right - Abell2495.5.13.0, a cluster normal SN la at 2 = 0.0796. Bottom loft 
- Abell399.3.14_0 is a normal cluster SN la at ^ = 0.0613. Bottom right - Abell85.6.08.0 is a normal cluster SN la at z = 0.0617. 
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Fig. 9. — The histogram of the velocity difference between the cluster redsliifts and our SNc la. Our hostless SNc la are shown in the 
filled histogram. Only those SN fa with At) < 8000 km s"'^ are shown. We choose a Av < 3000 km s^-"^ as our definition for a cluster SN 
la. Larger values of Av only add SNe beyond R200, beyond which MENeaCS becomes incomplete due to the Megacam field of view in any 
case. 
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Fig. 10. — Red sequence offsets of SN host galaxies in MENeaCS. The dashed lines indicate the median scatter in the red sequence over 
the enitre MENeaCS sample. Left SN host galaxy offsets for the entire MENeaCS cluster SN la sample. Ten out of the 19 SN la with hosts 
were consistent with the cluster red sequence, with bluer hosts becoming more prominent in the cluster outskirts. Middle SN host galaxies 
which were within 1 Mpc of the cluster center, of which 6 out of the 8 were consistent with the red sequence. Right SN host galaxies within 
^200 1 of which 9 out of 13 were consistent with the cluster red sequence. We use these SN numbers for determining out red sequence SN 
la rate in § [8] 
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Fig. 11. — Red sequence offsets of SN host galaxies in MENeaCS as a function of clustercentric distance in units of R200 (Left) and kpc 
(Right). The dashed lines indicate the median scatter in the red sequence over the enitre MENeaCS sample. Note that the most central 
host galaxies are all hosted by red sequence hosts, but the sample size is too small to say anything definitive. 




Fig. 12. — Sample supernova detection efficiencies (r7(m)) in the MENeaCS survey. Detection efficiencies for the central regions of our 
images sections (left) and corner image sections (right) are shown. Efficiencies for good seeing (black), acceptable seeing (blue), and bad 
seeing (green) are shown in each. See § [7] for details, along with an explanation for the detection efficiency fall off at faint and bright 
magnitudes. 
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Fig. 13. — A compilation of cluster SN la rates as a function of redshift over all galaxy types (left), and in early type/red sequence 
galaxies (right). Exact values, and their sources, are in Table|6] The star-poin t represents data meas ured in the current work. In the right 
panel, we also indicate the field elliptical SN la rate at 2 = 0, as determined bv lMannucci et al.l II2008I ) as the thatched box. The MENeaCS 
red sequence rate is consistent with the field elliptical rate at ^ = 0. 
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Fig. 14. — Recovered delay time distribution based on early type/red sequence SN la rates in clusters, taken from Table[6l and assuming 
that the constituent stellar population formed at z = 3 in a single burst. The star-point represents data measured in the current work. 
The solid DTD curve is the best-fitting power law, with en t" , and an s value of —1.62. The \a error region for the power law is 
— 1.08 < s < —2.16. We also plot the best-fitting and t~'^ power laws for illustrative purposes, as these are the expected power law 
forms for the double degenerate and double detonation scenarios, respectively. 



